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CHAPTER I 
INTRODUCTION 

It is recognized that there exist two modes of trans-— 
ferring heat between a surfuce and a flowing fluid. In the 
one mode, heat energy is transported by the fluid which is 
caused to move as a consequence of the heat energy it is 
transporting. The hotter and less dense portions displace 
overlying layers of colder and more dense portions and the 
fluid begins to flow, Natural convection occurs when heat 
is transferred by such a process. 

In the other mode, an external influence is introduced 
so that there is relative motion between the surface and the 
fluid whether heat is transferred or not. Forced convection 
occurs wnen hezt is transferred, 

When a heated pipe is supported in still air with the 
pine axis ina vertical position, heat is transferred from 
the pipe to the air by natural convection. The air moves 
emercd near the pipe well; the motion of the air is ina 
direction parallel to the pipe surface generators, Measure- 
ments by Koch “oe show that the diameter of the pipe has 
an effect on the hext transfer; the heat transfer increases 
as the pipe diameter is decreased. 

A conjecture on the effect of diameter upon the heat 


transfer in forced convection when tne fluid flows parallel 


to the surface generators has been mede by Jakob and Dow (2). 


5 References are correspondingly numbered in the Bibliography. 





Their theoretical reasoning, supported by measurements, 
shows that an increase in heat transfer may be expected as 
the diameter of the pipe is decreased, 

A distinction is required when heat is transferred by 
forced convection with the fluid flowing outside and parallel 
to the surface generators. One possibility is the case where 
the pipe is placed ina fluid stream bounded by the walls of 
a duct or other enclosure; the other is the case where such 
an enclosure does not exist. 

This distinction arises from the concept of boundary 
layers as envisaged by Prandtl (3). It is postulated that 
a boundary layer of fluid is formed on any surface in contact 
with the flowing fluid. The main characteristic of the 
boundary layer is that the velocity of the fluid within the 
boundary layer is less than that of the main stream outside 
the layer, becoming zero at the surface. If a pipe is placed 
inside a duct, a boundary layer will form on the walls of 
the duct and on the pipe surface. It is obvious thet the 
boundary layer forms only upon the pipe surface in the absence 
of such a duct or enclosure. 

When fluids flow inside pipes, it is known thet the 
boundary layer thickness eventually becomes equal to the 
pipe radius at a certain point, the point depending on spec-— 
ific flow characteristics, With this point fixed, the flow 
conditions are similar at places downstream from the point 
and all the fluid is in the boundary layer. At places up- 


stream from the point, portions of the fluid are inside the 








boundary luyer end the remuinder is outside the boundery 
layer; the relative proportions chenge from olace to place. 

Theory and experiment show that heat transfer for 
fluids flowing inside pipes is dependent upon the described 
boundary layer patterns, end it may be exvected that similar 
problems are presentec when the fluid is flowing outside the 
pipe. The presence of a duct or enclosure poses conditions 
that differ from the conditions in the absence of one, 

It is the ourpose of this work to determine the effect 
of curvature of a round cylindrical surface not surrounded 
by an enclosure, with the fluid flowing parallel to the 


surface reneretors, on the heat transfer by forced convection. 














CHAPTER I] 
THEORETIC:.L CONSIDERATIONS 
nen neat is transferred by convection, the process 


is mathenatically expressed as; 


og Pag 
a= = 1G, Orn. = == II-1 
as 9 dS 


This differential exoression shows thet the locel surface 
coefficient and the local temperature difference are taken 
ata particular point on the surface. 

The expression may be integrated when it is known how 


9 and h vary with S. This operation gives the total heat 


S 
q = Ic as II-2 
O : 


On the other nand, the total heat trensferred may be 


transferred as: 


written as: 
q 
q = h,@,S, or h, = --- II~% 
8.5 
{i] 
Substitution from Equation II-e into Equation II-3 


yields: 


1 
ho= -=- he ds TI—4 
o 


Inspection of Equation II-%3 shows that the choice of 
es will influence the value of hw for the oroduct hy @ 28 


constant for fixed values of q and S. 








C2 


1 
=~ 8 as II- 
: z) 


Substitution from Equation II-5 into Equation II-4 


Peis; 
S 
ho as 
QO 
hi oe tes 
8 das 
G 


The above developed relations are simple ones. The 
efiects of fluid properties, velocity of flow, and arrange- 
ment of surface are not readily apparent. Nor is the 
relationship among these effects apparent, Since Equation 
me. 1S used to define the local surface coefiricient, it is 
exoected tnat the coefficient is depencent upon these effects. 

The determination of these effects upon the coefficient 
is an enormous task. Hany variables are required to express 
the fluid proverties and tne surface arrangement. They may 
dnclude such properties as density, specific heat, viscosity, 
and therinal conductivity. The surface may be the interior 
mreeexterior of a pipe, or a ilat olate. There may be an 
adiabatic section of surface »receding. The flow direction 
may be parallel to the surface generators or at some angle 


tO them, 








The application of tne tneory of similarity to heat 
transfer problems by Nusselt (4) brougnt order from chaos, 
The procedures of dimensional analysis supplemented his 
work. In essence, these methods serve to assemble the 
many variables into dimensionless grouns and have the 
effect of reducing the number of variables to be studied, 

An excellent treatment of tnese methods may be found in 
Jakob (5). 

A oractical result of Nusselt's work may be illustrated 

by the equation: - 
Ny = f\Nper Nop» Q ) Tie 7 


that is, tne Nusselt number is expected to be a function of 


pet 


Reynolds nunber, Prandtl numver, and surface configuration. 
Other g¢rouos may be present. 

Tae anvdlication of dimensional analysis to the heat 
miiaster probleia requires that the form of the function of 
Equation II-/7 be a power function. This form is illustrated 
as: 


rT | n Dp 4M 
Nu Co (Nie) (Ny) p I1-é 


N 
where on n, p, and m ure determined by experiments. 
The inportance of Nusselt's work becomes more apvarent 
when it is realized that the values of Cs » nN, p, and m 
remain constant for almost all fluids when the flow config- 
uration remains the same. The major changes in these values 
occur when the flow configuration changes. 


In general, three flow configurations are recognized, 














They are termed laminar, transitional, and turbulent. 
Laminar flow along an immersed surface is characterized by 
the absence of a velocity comoonent normal to the surface. 
The direction of the velocity is always parallel to the 
surface. Investigation of the velocity shows that it has 
&@ value of zero at the surface and increases in accordance 
with the principles of viscous motion until the magnitude 
equals that of the main stream. It may be noted that a 
boundery layer has been indirectly described. The place in 
the fluid where the velocity of the fluid in the layer is 
practically equal in magnitude to that of the main stream 
describes a boundary layer thickness. 

There is no such orderly motion in turbulent flow. 
Except for a postulated thin lamninar sublayer at the surface 
of the body, the fluid supports velocity components ina 
Girection normal to the surface, and actual transport of 
material in this norial direction takes plece. Investiga- 
tion of the velocity component parallel to the surface 
shows that it increases with distance from the surface 
until it is practically that of the main stream. A boundary 
layer thickness is upain described. 

Transitional flow is the regime described when neither 
faminar nor turbulent flow occurs. This kind of flow is 
considered to exist in the region between laminar and turbu- 
ment flows. 

These concepts are more cleezrly illustrated in Figure 1. 


In Illustration (A), the velocity U,, and the length of the 








6 


Suriace oe ere such that a laiinur bouncery layer is 

fOTMmec on the Whole plate surface, If the velocity is 
incre:sea to ae a place L, occurs Wherein the boundary 
meer Orececing this piace 18 le&Miner, and a transition flow 
eccurs following tuis place, shown in Illustration (8). If 
moe velocity is increased further vom all three kinas of 
boundery liuyers may occur as snown in Illustration (C). 
eet retion (D) sLows how tue velocity may be distributed 
Mecnin tne bouncury layers, 

These concepts wre of reat importance in heat transfer. 
[the case of laminar flow, where &« velocity cornponent 
megial to the surface is lacking, heet energy is transferred 
Similcr to tne heet concuction process. In turbulent flow, 
Meere &. tr:nsoort of material normal to tne surface occurs, 
mee, energy is transferred along vito the material transport. 

im the orececings Giscussion of the velocity Cistri Duties 
Pipiiin the bouncdery layer, it was noted that the velocity of 
the fluic at the surface is zero. This means thet some 
particles of the Tluid heving the velocity of tne nain stream 
heve lost that velocity. The kinetic energy associated with 
the velocity has been aissipatead. This pnuenowenon is known 
@s aeroavnamic heating; tne texnerature of an adiabatic 
eurfece acquires a value gre:.ter then thet of the free stream. 
As tay be exnected, tunis effect requires attention in heat 
meenstfer problens. 

A descrintion hae been riven about the form:tion of 


Doundcery l@yers, and it nas been incietea iow the character 
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FIGURE 1. ILLUSTRATION OF BOUNDARY LAYERS 








se, 


of the layers may change. Another efrect muy arise, for 

the neat trensfer may or may not occur over the whole 
surface. Tne whole surface may be the heat transfer surface 
ea the heat transfer surface may be prececed by an adiabatic 
section or surface. These consicerations also require 
attention; for the bouncary layer at the place which marks 
the beginning of heat trensfer may ce of different thickness 
Mepchuracter, depending ujon the choice of the length of the 


auliabatic section oreceding the heat transfer section. 








PeevlOUS iuVesTiGuliIONs 
Tie description of the flow configuration is f¢iven most 
remiurkebly by the Reynolds number. Reynolds (7) end sub- 
sequent observers est«dlished the existence of usper and 
} 

lower critical Reynolds nutbers, (N_.). and (ih..)_, respect- 

: boes.C eae : 

' 

ively. For... =| ) a? she flow was found to be laminar; 


for = 


ae me TloW was found to be turbulent. Further 


igs 
' 

the velues of i) eG (Need are unique for all fluids 
Memeliy e:counterec. sy referring to Firure 1, this may be 
meLustrotec. is: 

V 1 4 t 
” ob obe ; 
( N- i See em, COs bane 


Wh ee 


Koreover: 


, L 
fia. = mone anotner constant. 
lols fs 


O 


The results of experrmental determinetion of tnese 


Peastants is ¢iven in Table i. 
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Tepes igi 
C27 TT LChie lO MS Vie POR FLAT FLATES 


\all duta recduceé to busis of total Length) 
’ 


Source ee) «. (Urea 
ie = 5 
Bo,.ert (1 ) a. 10 Fie eee 
Dodge «nd Thouwoson (&) -- 5 x 10; 
Paes (9) ~~ 4 7,6 x 102 
Ten Bosch (11) 1 x10 5 x 10° 
Pome oil, Hteovnsk Cal: wOose TIP 
4 
fakob and Dow (2) Sup ees. Re aes O52 


11 











A ives 


elias arrived at nis value by judging from mezsured 
velocity profiles. Ten Bosch, according to Seibert (72), 
La = 1x 10° may occur in the event of 
major cisturbance; that both critical Reynolds numbers may 


poates that (i 


be zero with a blunt leacing edge; and that the transitional 


< 
ops Nee = ex 10/ for sharpened 


4 


recion may exvengerrom | x 1¢ 
leacing edges and steadied flow, 

The selection of the characteristic length in defining 
Reynolds nuuwber is significant. Elias based his observation 


on the thermal length and judged (lip), = 3 x 10°. This has 


6 
tne effect of neglecting the boundary layer that has been 
formed on the preceding unheated starting length. The value 
from Hlias listed in Table I]@ias been corrected to include 
the starting; length; the starting length was ten centimeters 
and the thermal length was fifty centimeters. 

The value for the exponent p in Equation II-& is given 
Baeraole I]. 

TABLE II 
EXPONENT FOR PRaiiDTL NUP BER 


Laminar and Turbulent Flow 


Source Surface “xponent p 
KcAGams (6 Tubes and Plates O.4 - 0.3 
Colburn (14) Plates 1/3 
Pohlheusen (12) Plates Ws 
Latzko (13) Plates iL 


The values for p given by McAdams are the result after 
Sunmarizing the work of many investigators. He attributes 


the variation to viscosity effects erising from heating the 
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PPOTa iil One CaSemeaE@® Tromcoolimge in the other, and 
recomuends & value of 1/3 for fluias of low viscosity or 
if viscosity effects «wre considered by eddition of another 
term to #quation II-c. Colburn's value is an assumed value 
erising from an effective correlation of heat transfer data 
anc comparison with fluid friction. Pohlhausen's value is 
theoretically derived for laminar flow. Latzko's value is 
tueoreticully agerivea for turbulent flow over the whole plate 
Puriacce. It may be noted that if the Prandtl number is near 
unity, then the effect of the exnonent is small. 
Values for n in Equation II-& are given in Table III. 
Peet III 
@XPOUSNT FOR R£EYUOLDS NUMBER 
pource surface Exponent 


Laninar Flow 


Colburn (14) Plate 0.5 
Pohlheusen (12) Plate 0.5 
Jakob and Dow (2) Cylinder Ons 
Prener anc Norris (18) Cone end Cylinder O25 


Turbulent Flow 


Colburn (14+) Plate Ors 
Juerges (19) Plate 0.775 
klias (y) Plate Greater than 0.4 
Latzko (17) Plate CZs 
Griffiths anc Awberry (20) Cylinder Ta 
Jukob and Dow (2) Cylinder 0.8 
Meener and iorris (18) Cone wnd Cylinder 0.8 


mererrin: to Table Iii, Coiburn's values result froma 
Sumiary of various investiputors work. Polllé~usen's value 
eemes frou tueoretical consiaeretions. The values revorted 


‘ 
by Hlias and Griffitus and Awberry are orovably cue to a 
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li.ck of fully Geveloped turbulence. The data of the other 
observers listed are reesonably well correlated by the 
Penorved value for n. 

The selection of the characteristic length for defining 
Reynolas nuuber is not significant with respect to the value 
of n. The transformation of the characteristic length from 
thermal length to total length and vice versa, for example, 
igs acconplisneda by &@ length ratio wnich does not affect the 
Mmereuc Of the exponent. To illustrate, let 


Hh 
pd a NY 
‘ = f e e r e 
Bey 6 (Noe dy 


i wa TeV.) 
and Gi), =o Ga on. 


fe) 


Q 

Oo = 
> 
I< 
<I 
O° foe 
<i 
nN 
= 


r ue X\nf/L a 
or (Hwy, a0 ae SG ee 
. Vv, x 


Therefore, the exponent for Reynolcs number is the same 
Mmerespective of the cheracteristic length used in definition 
as long as the selection iS consistent for any given case. 

Investigation of the function 9” is a relatively recent 
occurrence. The first conclusive evidence of considering a 
starting length function is offered by Jakob and Dow (2). 
Using the mean surface coefficient, and the total length in 
both Nusselt and Reynolcs numbers, they arrive at tne 


exoression: 


ae lo aie oe Ofion = (ay) = 02606 
Thea 
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for turbulent end leainer flow. This result was determined 
exoserinentally by neat transfer measurements with unheated 
Sterting leniths. also, they offer a result wnen the 
.usselt iumaber 18 com uted witn the local surface coefficient, 
Beecise: by cifferentiation. Witi certain restrictions, tney 
meport: 

B= 0.5 + 0.2(8/L) ~ 0.78(6/L)*"'? + 1.18(s/L)-° 

IiI-2 
Bnoticr exoeris@ntal Tesult Was recently oiter ome, 

faisel und Snerwood (16). Using the local eurface coeffic- 
Ment, «end total length in both nusselt and Reynolds numbers, 


miey report: 


id 


g? = fr - (9/1) ®|> 9-44, 0,us £ (s/t) £ 0.95 
Ti9 


for turoulent flow. Tiuis result comes from investigating 
Wass transfer and the starting lengtn may be consicered to 
me cClabatic. 

Theoretical results have been revorted. Using the 
local surface coefficient, anc total length in both Nusselt 
and Reynolcs numbers, Eckert (10) and Rubesin (15) each 
mepOrTt: 

gu = E 2 (6/1)2/"4 eee, for All ee 
III-4 
mer laminar flow. 

Ph thneoretical result for turoulent Plow is cue to 
Rubecin (15). With local coefficients und characteristic 


Menethis as before, ue gives: 
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gr = E - (6/1)39/40 | _ i ee forall (s/L). 
III-5 

Comparing these equations, it is seen that the express- 
ions in Equations III-3, III-4, III-5, tend to becone 
infinite as (s/L) becomes unity, whereas the expression in 
Equation III-c aporoaches the finite value 1.40. Concerning 
Equation III-2, Jakob and Dow (2) state: “However, since 
differentiation of an empirical equation is a celicate 
matter, the accuracy of the result should not be overestim - 
ated.". 


HE 
The expression oe of Equation JI-%8, may be 


Nop)? 
given by another constant om, which assumes different values 
according to the fluid properties. Values of C. fOr aura 
atmospheric pressures and about 68°F are given in Table IV, 
according to Jakob and Dow. The values in the table are 
corrected for the ratio (s/L) according to Equation III-1l. 
TABLE IV 
VALUES OF C, FOR AIR (FROM JAKOB AND DOW) 
Atmosoheric pressure and about 68°F. 


Source surface 8 


Laminar flow, n = 0.50 


Pohlnausen Plane Ooty 
Fage and Falkner Plane Oi 28) 
Jakob and Dow Cylinder Oe 570 


Turbulent flow, n = 0.80 


Latzko Plane OO” Oe54 
Juerges Plane O07, Ogee 
Colburn Piene 0.0320 
Seibert Plane 0.0349 


Jakob and Dow Cylinder 0.0280 
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Reeults recently obtained by Berman (21) may be added 
to this table. ne reports o = 0.3/ for laminar flow. 
This result comes from mass transfer measurements from 
Cylinders of ice having hemispherical nosetips and with 
cylinder ciwmeters from e.ee inches to 3.79 inches. 

The efrtects of aerodynamic heating on heat transfer 
have been studied and a summary of results is given by 
Johnson and Rubesin (17). The sum:ary shows how heat trans- 
fer results obtained at low fluid velocities, wherein the 
effect of aerodynamic heating may be neglected, can be ap- 
olied to measurements where tne flow velocity is large and 
may be greater than the speed of sound. 

A sumuary of the foregoing brief review of previous 
investigations is presented in Table V. 

TABLE V 


SUeutAKY OF PR2zAVIOUS INVeS5i IGATIONS 


guantity AD )roximute Value 
(tre) 5 x 10 
“Re/e oie ae. 
p 1/3 
n, laminar O50 
n, turbulent 0.80 
Gy Discussed below 
Co tf uN 


r W iT 


Altnou;n tne function @ nay avoear to have inconclusive 
TOrms, it is noted that Laisel and Sherwood reported that 
the Jakob and Dow fom fave a cood correlation for their 
data ulso. Chis oroblem will be discussed in wore detail 


later. 
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When one of the groups on the right side of Equation 
II-& is held constant, the effect of that group is absorbed, 
civing a new equation constant. This was illustrated in 
the presentation of Table V in which the expression C (1 2ee 
was replaced by Co. A similar effect may be expected with 
respect to a curvature factor. 

The Jakob and Dow result for laminar flow agrees with 
a theoretical result for a plane. However, the theoretical 
result for a plane has an intermediate value between exper-— 
imental results for a plane and a cylinder. Therefore, the 
effect of curvature is not too clear, although the tneoret- 
ical reasoning of Jakoo and Dow indicctes that results for 
a plane may apply to results for a cylinder of 1.%3 inches 
erameter, 

For turbulent flow, the Jakob and Dow result is between 
a theoretical one and experimental ones. If their theoret-— 
ical reasoning regarding curvature is correct, then there is 
strong indication that the theoretical value for a plane is 
nearly correct. If their theoretical reasoning is incorrect, 
then comoarison of the experimental values indicates that 
the effect of curvature may be to reduce the heat transfer 
4s curvature increases. Therefore, tne effect of curvature 
is in doubt for both turbulent flow and laminar flow. 

The above considerations indicate a plan for exoeri- 
mental attack; 
(a) the curvature must be made great enough so that a com- 


oarison with the flat plate and a cylinder 1.3 inches diameter 





c) varying thermal lengths should be used to determine at 


a - 
" NO wrol) 
AD wo Ad 


the function ~”° is valid for variable thermal length | 


48 variable starting length. 














CHAPTER IV 
APPsRATUS AND INSTRUMENTATION 

It is desirable to express heat transfer data by a 
relation as iquation II-é, for the equation is based upon 
the powerful tool supplied by the theory of similarity. 

The application of such an equation to the flow of fluids 

in pipes, for example, shows the effectiveness of the theory 
over a very wide range of variables. Consequently, Equation 
II-S is to be used to express the results of this work, 

Proceeding term by term, it is seen the Nusselt number 
requires measurement of a surface coefficient of heat trans- 
fer, @ characteristic length, and the thermal concuctivity 
of the fluid. Reynolds number requires the relative veloc- 
ity between fluid and surface, a characteristic length, and 
the kinematic viscosity of the fluid. Prandtl number re- 
quires the specific heat, viscosity, and the thermal conduct- 
ivity of the fluid. The function @ requires two character- 
Peete lengths of the surface. 

The measurement of the surface coefficient of heat trans- 
meeemay not be made directly. Referring to either definition, 
Equation II-1 or Equation II-%, it is seen that a knowledge 
is required of the heat transferred, the temperature differ- 
ence between fluid and surface, and the heat transfer surface 
area, 

Fluid properties as thermal conductivity, specific heat, 
viscosity, kinematic viscosity, density, and others are 


required. A knowledge of any two properties, such as 
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temperature and pressure, serves to determine all others. 

Rernoulli's Theorem affords a simole method to deter- 
mine the fluid velocity. <A knowledge is required of the 
mead oressure, specific weight and the local pravitagiona) 
constent. 

The heat transfer apparatus is illustrated in Figure 
e. Heat is transferred from the long slender element to an 
air stream flowing parallel to the long axis. A hemispher- 
ical nosetio is observed at tne upstream end. The down- 
stream end is clamped to a vertical stanchion which is 
mounted on a wood base. The stanchion is braced by guy 
wires attzched to the base. 

Vibration insulation of the base from the floor of the 
laboratory is accomplishea by the use of shear-type rubdber 
Shock mounts. A vibration ebsorber is provided by the 
weight and rod fastened to the stanchion. These vibration 
controls are illustrated in Figure 3. 

The heat transfer element consists of an outer tube, 
an inner tube, an electric heater, and selected nosepieces. 
The outer surface of the outer tube vrovides the heat trans- 
fer surface. The inner tube, which slidaes into the outer 
tube and surrounds the heater, is used for fastening the 
nosepvieces to the outer tube. 

The outer tube was fabricated from brass tubing, seventy 
per cent copper and thirty per cent zinc. Three lengthwise 
holes, one-sixteenth of an inch square, located at 120° in- 


tervals, were provided for insertion of thermocouples. 
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Fabrication of the outer tube is illustrated in Figure 
4, The tube was annealed by heating at 750°F for one-half 
hour, prior to machining. The machining process consisted 
of cutting grooves one~-sixteenth of an inch wide by three- 
thirty seconds of an inch deep followed by cutting second 
grooves, centered over the first, one-eighth of an inch 
wide by one-thirty second of an inch deep. The upper groove 
Surfaces and strips of brass rod, one-eighth of an inch 
square, were tinned with soft solder - sixty per cent lead 
and forty per cent tin. The strips were placed over the 
grooves and clamped to the tube at two-inch intervals. The 
sweating operation was accomolished by using an electric 
heater inserted into the bore of the tube. After cooling, 
the protruding portions of the brass strips were machined 
away. The surface was highly polished, bright chromium 
plated, and finally repolishec. 

The interiors of the square holes were cleaned and 
polished. This was done by sliding and rotating steel drill 
rods, coated with valve grinding compound, in the holes. 
Another piece of drill rod was used as a gauge and the holes 
were judged to be clean when the gauge rod could slip easily 
through the holes, 

It is realized that the continuity of the sweated joints 
is of great importance, and that a thermal bond must exist 
along the whole joint. It is reasonable to assume that a 
thermal bond was obtained for the reasons following: 


(a) It was necessary to perform the sweating operation many 
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(a) First Groove (bd) Second Groove 
(c) Tinning (d) Sweating 





(e) Completed 


FIGURE 4, FABRICATION OF HEAT TRANSFER SURFACE 
(SCALE: 2"=1") 
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times before it was found how to prevent the solder from 
filling the holes, Each time the parts were separated, the 
Surface showed that a good bond had existed. 

(b) A viece of finished tubing, about one inch long, was 
cut from one end. The cut was made at the middle of the 
first clamp interval. Mechanical destruction of this piece 
Showed complete bonding had occurred. 

The outside diameter of the finished outer tube was 
found to be 0.624 in., plus or minus 0.002 in. The inside 
diameter was measured to be 0.%54 in. The tube was thirty- 
five inches long. A surface plate was used to determine 
longitudinal eccentricity; this was found to have a maximum 
of one-sixty fourth of an inch located eleven inches from 
the end subsequently used as the downstream end. 

The inner tube was of brass, seventy per cent copper 
ana thirty per cent zinc, und was not annealed. This tube 
had an outside diameter of 0.%76 in., and an inside diameter 
of 0.326 in. One end was closed with a brass plug, one- 
Cignth of an inch thick, silver soldered to the tube. The 
olug was drilled and tapped for a Number 10-3e machine screw, 

At the plugged end, the outside diameter of the tube 
was reduced to 0.246 in. for a distance of one-fourth of an 
inch by cutting away the metal. This device serves to reduce 
heat conduction to the end plug. 

Tne plugged end of the inner tube was designed to extend 
one-fourth of an inch beyond the upstream end of the outer 


tube. The other enc of the inner tube extended beyond the 
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downstream end of the outer tube and wus flared to fit a 
compression tube-fitting assembly. This assembly served as 
a device for tightening nosepieces to the outer tube. A 
set screw was provided in the assembly to anchor the elect- 
mac heater. 

The electric heater was a General Electric Company 
tubular Calrod heater (22), rated at 1500 watts, 230 volts. 
This element consists of an outer heat-resisting steel 
sheath surrounding the heating element wire which is embed- 
ded in magnesium oxide insulation, 

Initially, the outer sheath was forty-eight inches 
long between terminals and was insulated from them with mica 
washers. The active heating length was located within the 
middle thirty-six inches, leaving six inches at each end un- 
neated, 

wodification of the heater was accomplished by cutting 
ten inches from an end. This cut exposed the heating element 
cross section. Here, the outer sheath was cut and filed 
longitudinally to form four gores so that a closed hemispher- 
ical tip could be shaped at the cut. The wire of the intern- 
al heating coil was placed into these cuts while forming was 
in progress. The wire was firmly pinched and held by elastic 
forces to the outer metal sheath. Care was taken during the 
shaping process to eee loss of insulation, 

This construction served to make a heating element with 
one terminal as originally provided and with any place on 


the outer sheath to be used for location of the other 





or to comolete tne | 


As woGified, u terminal lug was solaered to the 


adjacent to the other terminal, 

The vopbion of the Calrod heater cut off was usedeso 
establish pertinent dimensions for the heeter. These were — 

s follovs: 

Cutside Giumeter of outer metal cheath 0.7316 

Insice diumeter of outer metal sheatn Ove oe 

Pectingwcolmeyvere cimmeter 0.020 

Unhe: ted length 6.25 

From touis inrornation and by measurements, the modif- 


Meo heater was thirty-elgnht inches long, The active heat— 





ing length was tuirty-one and three-fourths inches long, 
Derinning at tne formed end. 

Five noseoieces were made of Super maple aqowel ‘Stee 
fe hed nemispherical nosetinos. They were carefully cut 
and sanded so thi.t the finul cismeter was the same as that 
Meetiice outer tube. 

A hole, twenty-three sixty-fourths of an inch in 
Seameter oy one wnd one=jhulf inches deen, wes crilled into 
mie beses of four nosesieces. another hole, one-eirhth of 
Meeernch in @isameter, wes drilled along & nose viece aadmeter 
Bee i::ch from the bese. aA glued maple pin was driven into 
Peis hole aac throujin the eye of a linking device inside 
tue nosepiece. Tne pin was Cressec to size and smoothed to 


mac nose xyiece curface. 
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A lip cut into the nosepiece bi.se served as a locating 
Gevice, The lip enguped into a rece@e in the outer tute, 

The linking Gevice wus mude of brz«ss vire, one~sixteenth 
Of an inch in dicmeter. An eye was formed at one end and 
the other end was silver soldered to a brass stud. The stud 
Wee threadec to fit the plug in the inner tube, 

the four nosepieces constructed this way were five and 
oue-fourth inches, twelve and one-sixteenth inches, eighteen 
and one-sixteenth inches, «nd twenty Tour inches long, res- 
pectively, measured from tie joint with the outer tube to 
Mere tanpent plane ut the nosetip. 

ene fiith nosésiece, five-sixtecnitis of an inch Loney 
regquirea diiierent construction. The hole in the base was 
One-—cighth of an inch deep end a one-eighth of an inch hole 
mee counteroored to a centh of one-eighth of an inch. The 
Sinking; device Was a orass stud with a Number 5¢e hole 
mee lied along @ diameter at one end. <A corresponcing hole 
me crillec along @ nosépiece Gailameter and the linking de= 
Vice assembied with a brass Wire. 

The aetalls of tnese asseublies are illustrated in 
meeure 5. 

Uoon assembly of the nosepieces to the outer tube, the 
Meting surfaces were couted with Duco tiousehold Cement. 
After asseuioly, when the cement hed dried, any cement that 
mea oozed from the joint was reuoved. The nosepiece and a 
Mertion of the outer tuoe Were fiven two corxts of lacquer 


and tnen smwovtned with fine steel wool and rubovoed witn Duco 
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Rubbing Compound until a smooth joint was obtuined and no 
mecaquer reimiuined on the outer tube, 

wach joint was tested for smoothness. This vias done 
by placing the assembly in V-blocks upon a surface plate. 

A macuinist's dial-indicator was used to test the joint. 

My joint imperfections were treuted with lacquer and re- 
polished until the dial indicator showed zero deflection 

Bee the joint. 

Beiore assecmoling the electric neetver Go titewemocs 
tube-nosepiece assembly, a mark was mude on tne heuter sur- 
face. This mark served to set the insertion depth so that 
the heater tio woulda be in the olane of the outer tube- 
nosepiece joint. The heuter was inserted to the estiblished 
mark and then locked in place with a set screw, 

The heat transfer element was electrically insulated 
from the stanchion. <A smoothly bored one-inch IPS pipe tee 
at the too of the stanchion was split horizontally and each 
helf was lined with insulating paper. A finned spacing de- 
mree wes used as @ filler piece between the neat transrer 
Meement ana the pipe tee. The finned Somcine device alee 
served as a cooling: elewent. Hose clanps were used to hold 
the assembly. 

Direct current was used as the source of heat energy. 
The current wus providea by a motor-generator set equipved 
with a General Klectric Comwoany wlternating current induct- 
don motor, rated et fifteen norsepower at 1cO00 revolutions 


per minute. The motor was connected to an Electro Uynamic 











Comoany alrect current generator rated at 0-110 volts, 115 
amperes, at 1200 revolutions per minute, and to a Leece-— 
Neville Cowoany twenty-four volt aircraft—-tyoe generator, 
This latter twchine was used as an exciter. 

An aircrzft-typve voltage regulator was modified so 


ve 


gaat the voltage outout of the main @enerator could be con-= 
Mmeottec over the full rated voltae range, 

The r@latively hich mechanical inertia of the motor- 
@eneravor set in conjunction With the modiried voltage 
Meeulator fuve ood regulation. The voltage variation cid 
me Exceed one-hulf of one per cent for the current magni- 
tudes used, 

The heater power consumption was obti.ined by measuring 
the voltare crop between the hee«ter terminals und by measure 
mae the electric current. The voltage droo was measured by 
mastalling a moaified Leeds and Northrup Cornpany volt box 
maepare2llel with the heater. The current was measured by 
using a Weston Instrument Company shunt in series with the 
meer, 0th measuring elements were connected to a Leeds 
ema viorturuo Company oortable orecision notentiometer 
through a rotary switcn. 

A blower and auct system served to provide an air 
stream. ‘The blower was @ Power xngineering Company Tyne HS 
Meenine, driven through a Reeves Vari-sSpeed drive, by a 
mive horseoower induction motor, 

The straight duct system wis connected to the blower 


by means of & transition piece and termineted with a nozzle 








SLoWer Giccherme flance anw the nozzle mouin. 


was nine end seven-eifnhtns inches in ciameter. 
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Control or air velocity was ootwined with the Reeves 
drive and ea set of removable cGumoers at the blower inlet. 
‘The blower inlet, equinsed with a@ osychrometer, was in the 
laboratory roon. 
Tame measurement of velocity hea 
Pere impact tune nlaced into tne Cuct three 
Two Ellison Draft Cages (2%) were used for 


andicating the velocity heaa. The low-renge frage was used 





mor Velocity creacs below. @.5 in. of Water. The scale’ més 
Mmm inches long gra@@etec to 0:005 in. of mater. The hie 
Menge cepe, to fiwe incnes oI water, Was thirty inches long 


meecuated to 0.U1l in. o: Weter. 


} Measureéiiient of the @ir stream temM¢crature and estime-— 
} Mia of tne heat tranefer surface temperature Was accomp— 


bismed with seven thermocouoles. The thermocouples were of 
BMeecs anc Northrup Comgany coonver and constantan wires, 
Humber 30 Brown and Skarpe Gauge, 1736 Calivration. 


O ingicéte alr temper-— 


ct 


rour there@eouvles were used 
meure. Tnese thermocouples Were enazellec, single cotton 
mere reca wires, connected in series. Tite hot junctions*Were 
located at four equidistant places arounc the nozzle circun- 
Mmerence, Trey extended @m§to tne air streemh one inch away 


from the wall and were varallel to tne uir stream a aistance 








Se. one and one-half inches pointing upstream. The cold 
junctions were inserted into glass tubing wells imsersed 
mea mixture of melting ice und Water in Dewar tines. 

The conbinea thermo-electromotive torce of the four 
thermocounles was imeusured with a Leeds and Northrup Company 
Speedomix Indicator-Recorder. The indicator scule was 
eavided into tenths of a millivolt for a range of six midi 
moots. A vernier@was fabriceted and mounted on the poinzger 
mo perniit részaing@ to G.O1l miilivolts., 

mn Sweitcnine arrangcenent orovided a meéane to chee 
epeecdouex operation witn the sortuble precision sotentio-— 
moeuer. 

Three indepenaent thermocouples were used for estima-— 
mane surface temperature. These thermocouples were enam— 
elled, double wrap glass Tiber covered wires. They were 
fisted witn a onejw-hulf inch pitch and then dipped in iacq= 
uer and dried several times. Paint marks were adaed at 
intervals of two inches beginning at the hot junction, and 
Mae Ginding ena drying process repeated. The dipoing oro- 
cess served to stiffen the wires for eusy insertion into 
the heat transfer element holes, The cold junctions of 
these thermocouples were placed in an ice-water mixture as 
oreviously described. 

A ten-point Leeds and Northrup Company rotary switch 
wes used for selectively connecting the electromotive forces 
from the volt box, shunt, Speedomax, and the three thermo- 


couples 2avove to the pvortable precision potentiometer. This 
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instrument was provided with two scales. Electric poten- 
tials of sixteen millivolts and less could be read to 0.01 
nv. and estimated to 0.001 mv. -Potentials greater than 
sixteen millivolts could be read to 0.1 mv. and estimated 
to 0.01 mv. 

A mercury barometer was used to determine atmospheric 
pressures, 

An impact tube probe device was used to locate the 
heat transfer element in the air stream. The device, illus-— 
trated in Figure 6, was made by soft soldering two brass 
tubes to a split sleeve. The tubes were located one-half 
of an inch and one inch, respectively, from the surface and 
Berallel to it, pointing upstream. <A draft gage served to 
indicate velocity head. The sleeve was made to slip over 
the element and was provided with a clamp. 

The heat transfer element was judged to be centrally 
located and parallel to the air stream when rotation of the 
Gevice around the element showed no change in the velocity 
head. The device proved to be highly sensitive when located 
near the supporting stanchion as shown in Figure 6. The 
device was removed when making heat transfer measurements, 

A system for measuring vibrations of the heat transfer 
element was orovided. This system consisted of a motion 
picture projector equipoed with speed control, erected to 
cast a shadow of the heat transfer element tip upon a ground 
glass screen. The system provided a magnification factor 


of 2.9. Amplitudes of vibration were measured from the 
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screen and frequency was measured with a tachometer driven 
by the projector. 

The apparatus enabled experimentation within the 
following limits: 

Heat transfer element diameter 0.624 in, 


Nosepiece lengths O.31, 5.25, 12.06, 18.06, and 24.00 
inches 


Thermal length to 30 in. 
Ratio, x/L 0.14 to 0.98 


Air stream velocity 9.9 to 123.6 ft sect ; 
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CHAPTER V 
PRELIMINARY TESTS AND INVESTIGATIONS 
The method of measuring power consumption has been 
described, and it was noted that a modified volt box and 
shunt were used in conjunction with a potentiometer. The 
volt box was modified and tested as described in Appendix 


I, with test results as: 


Range, volts Multiplier, volts per millivolt 
O- 15 0.41992 
QO = 150 1.992 
O = 300 3,984 


The shunt was calibrated against a standard resistance 


as described in Apvendix II, with test results as follows: 


Range, amperes sUltiplier, amperes per millivolt 


The thermocouples were calibrated as described in 
Appendix III, and were found to be within 0.3 F of the cal- 
culated temperature of condensing steam at prevailing atmos-— 
pheric pressure. The three used in the heat transfer ele- 
ment differed from the others by 0.e¢ F. The thermocouples 
were calibrated in place after all wiring connections were 
made so that any effects of switches and lead wires might 
meminciucded. 

The electrical resistance of the modified heating 
element was measured with a Wheatstone bridge with results 


as follows: 
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Part Resistance, ohms, at 90°F, 
Heating coil 32, 34 
Outer sheath 0.060 


These data show that the error made by considering all the 
heat as coming from the coil may be negligible. Also, since 
the heating element was not insulated from the inner and 
outer tubes, the electric resistance attributed to the outer 
sheath may be reduced to some smaller value because of poss- 
ible parallel paths for the electric current. 

The impact tube located in the duct was compared with 
a pitot tube of Prandtl's design (24). During testing, the 
pitot tube was located twelve inches downstream from the 
nozzle mouth in the center of the stream. Assuming that 
the velocity coefficient for the pitot tube is unity, it 
was found that the velocity head obtained with the impact 
tube multiplied by 1.03 would equal the velocity head ob- 
tained with the pitot tube, All subsequent measurements for 
stream velocity were made with the impact tube in the duct 
and corrected accordingly. 

An exploration of the jet without the heat transfer 
element in the stream was made with the pitot tube. Down- 
stream stations were selected at six inch intervals begin- 
ning one inch from the nozzle mouth. Traverses were made 
at intervals af one inch horizontally and vertically at 
these stations. This »rocedure was followed in detail for 
a@nominal jet velocity of 10/.8 ft per sec. Explorations 


in less detail were made at other velocities. The results 








HO 


of these explorations for a central core of four inches in 


diameter are: 


Distance from Nozzle mouth, Loss of velocity, 
inches Der Cen, 
37 O 
4 1 
4g 5 


The exploration was repeated with the heat transfer 
element in place and with the tip located five inches down- 
Stream fron the nozzle mouth. The results for a central 


core four inches in dlameter are: 


Distance from Nozzle mouth, Loss of velocity, 
inches per cent 
af 5 
4g 10 


As a consequence of these data, all heat transfer measure- 
ments were limited to the portion of the jet included 
between five inches from the nozzle mouth to forty-three 
inches from the nozzle mouth. 

The Calrod heater was tested for uniformity of heat 
generation prior to modification. An estimate of the un- 
meormity was obtained by inference from surface temperature 
measurements. Tne method is based on the suoposition that 
the surface coefficient of heat transfer in free convection 
may be assumed to be constant for small changes of surface 
memoerature. Therefore, Variations in surface temperature 
indicate variations in heat generation as expressed by 
Equation II-l. 


These tests are described in Appendix IV and the 
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results are illustrated in Figure 7. It is noted that the 
Calrod heater was hotter at one end. It is assumed that 
the test results are expressed reasonably well by a linear 
relation with respect to length. The equations of such re- 


lations, obtained by the method of least squares, are: 


QO, = 307.1 - 0.43 z, F, at 125 watts input. 
(Appendix IV, Part A). 
63 = 160.4 -~ 0.19 z, F, at 41.2 watts inout. 


(Apoendix IV, Part B). 
These equations were derived for data within the middle 
twenty-four inches of the Calrod heater. 
It was noted in Chaoter III that heat transfer results 
Tor forced convection with turbulent flow between a flat 


plate and an wir stream may be exoressed as: 


a. 0.60 
(Mp = So NRe yp 


7 : ue 0.850 
or hie ae (k/1L) oe 
' 0. &0 0.20 
=O (Ne) a) (1/1) 


It is seen, therefore, that the surface coefficient becomes 
smaller as L becomes larger, other conditions remaining un- 
changed. Hence, if a plate is constructed so that the heat 
transferred per unit area is constant over the whole area, 
tnen the temperature difference between the plate and en- 
vironment must increase as L increases, 

As a consequence of the above discussion, consideration 
of the temperature difference of the Calrod heater as tested 


and the plate as predicted resulted in the conclusion that 
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the modified end of the heater should be the hotter end. 
A partial compensation for temperature difference in order 
to approach an isothermal heat transfer surface was effected 


when the heater was cut accordingly. 
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FIGURE 7. CALROD HEATER TEST DATA. SOLID LINES ARE LEAST 
SQUARES LINES FOR 12 - 36 INCHES, INCLUSIVE 





CHAPTER VI 
HEXPZRIMENTATION AND DATA 
The experiments were classified into seven series of 


runs as follows: 


series Nosepiece Length, Angle to Stream, Vibration 
Designation inches Geerece controls 
A Cel O Yes 
B Diao 0 Yes 
C 12.00 0 Yes 
D 18.06 0 Yes 
E 24,00 O Yes 
X 5.25 3 Yes 
V 5.) 0 No 
A run of any series consisted of making a temperature 


traverse along the heat transfer element at intervals of 
two inches after a steady state was judged to prevail. 
Power input and air stream velocity were the controlled 
variables for a given series. Variables as hunidity, air 
temperature, and barometric pressure were uncontrolled 
variables. The steady state was judged to prevail when the 
temperature difference between the air and the nosepiece 
Heant remained steady. 

A run was not started or was discarded if the Speedomax 
record of air temperature showed instantaneous variations 
exceeding 0.5°F or if the average air stream temperature 
risé was in excess of 1°F in ten minutes. 

The temperature traverse was made by moving the thermo- 
couples from place to place along the heat transfer element 
length beginning at the nosepiece joint. Preliminary tests 


revealed that the thermocouples reacned relative equilibriwun 


nn 
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within one minute after being moved. To provide an ample 
margin, data for the temperature traverses and correspond- 
ing air temperature were taken at intervals of three min- 
utes. 

It may be demonstrated that the surface temperature, 
the heat conducted to the nosepiece, and the accuracy of 
the determination of temperature difference increases with 
power input, whereas the accuracy of the surface tempera- 
ture measurenent decreases with power input. It was shown 
previously that if the unit surface heat generation is con- 
stant, a rise in surface temperature with length may be 
expected. On the other hand, a thermocouple does not yield 
reliable results unless the wires are in an isothermal zone 
for a distance depending upon the wire size and temperature 
gradient. As a compromise, therefore, the power input for 
a run was adjusted so that the maximum temperature gradient 
along the heat transfer element in the region of the nose- 
piece joint was in the neighborhood of one degree Fahrenheit 
mer inch. 

In all except Series V runs, vibration control was ob- 
tained by changing the absorber weight or length. Ampli- 
tudes could be limited to less than one-sixteenth of an inch 
in most cases and were about one thirty-second of an inch 
generally. The vibrations encountered were of sporadic 
nature and probably caused by nearby vehicular traffic. In 
the Series V runs, the shock mounts and vibration absorber 


were removec and the base was clamped to the laboratory 
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PLOoOr. 

Readings for power input determination were m:de four 
times per run and other readings were made once per run ex- 
ceoting surface and wir temperature measurements. 

The basic data are given in Tables VI to X, inclusive, 
and represent ninety-two runs. A total of ninety-five runs 
were made ana three runs were rejected because of failure 
to take complete data. 

The data of Table VI, Imoressed Heater Voltage, are de- 
rived from the 0 - 150 volts range of tne volt box. . 

ire Gate in Table VII, Electric Heating Current, are 
Gerived from the 0 —- 5 ampere range of the shunt and in- 
cludes the current passing through the volt box. The data 
of Aooendix I show that the maximum current passing through 
the volt box is 0.005 amperes with rated voltage impressed 
on the terminals used here. Examination of Table VI shows 
that the maximwn impressed voltage was 43.6 volts. The 
corresponding current is (43.6 x 0.005)/150 equal to 0.0015 
amperes. By Table VII, the total current for this voltage 
is 1.3385 amperes of which 0.0015 amperes is one-tenth of 
one per cent, aporoximately. A computation with the minimum 
values yields a Similar result. Tnerefore, no correction 
was made and the data in Table VII are used directly. 

A check between correSponding values of Tables VI and 
VII is provided by use of Ohm's Law and the measured elect- 
ric resistance of the modified heater. 


The data in Table VIII, Initial Air Stream Velocity, 
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are the result from considering the observed impact tube 
velocity head corrected by the multiplier previously des- 
cribed, corrected barometer observations, and the air stream 
temperature. To facilitate computations, the conversion 
tables of the Ellison Draft Gage Company (2%) were used. 


The tables are based on the formula: 


The maximum relative humidity at the blower inlet was 
thirty-five per cent, hence the air was assumed to be dry 
and no correction for humidity was made in the velocity de- 
termination. 

The temperature data in Tables IX and X are uncorrected 
for the steam point deviation found in Appendix III. The 
agreement between the thermocouples used for air temperature 
measurement and the computed steam point temperature is such 
that correction is not required. 

Comparison of the air temperature thermocouples with 
those used in the heat transfer element shows a deviation 
of 0.2 F at the steam point. Examination of Table X shows 
a maximum temperature difference of thirty-one degrees Fah- 
renheit, approximately. If it is assumed tnat the deviation 
of the theraocouples is a linear function of the temperature 
above the ice point, then the probable deviation between the 
thermocouples may be a maximum of (31 x 0.2)/160 equal to 
0.034 F, and is considered to be negligible. 


Of the data in Table X, Temperature Distribution Data, 
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it is assumed that the measured wall temperature is ata 
place corresponding to half the outer tube wall thickness. 
No correction is made for conduction of heat along the 
thermocouple wires because of the small wire size and the 
relatively small temperature gradient. The computed wall 
temperature is the average of the observations from the 
three traverse thermocouples. 

The Leeds and Northrup Company Temperature Conversion 
Tables, Standard 31031, were used. For purposes of inter- 
DOlation, the ratio of temperature difference above the ice 
point to the corresponding potential difference versus po- 
tential difference was plotted. A smooth curve was drawn 
through the plotted points and a new conversion table pre- 
parea for immediate purvoses. A cross check between the 
standard established the acceptability of the new table, 

Some of the data of Table X are shown graphically in 
Figures &, 9, 10, 11, and ie. Runs having nearly the same 
velocities were selected for the figures so that a visual 
comparison may be made. 


Mibration data, gathered from the Serres V runs are. 


Run NOS. Amplitude, Revolutions 
inches ner minute 
2 aly aie HUG 
7,45 0 7 
Saves 1/8 4 
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Comparative data from the Series B runs are: 


Run Nos, Amplitude, Revolutions 
inches per minute 
1,2,3,4,5 iy ate 520 
PY acy, 1/32 520 
HO, 11, 12,2 5,14, 15 1/6 520 


The data for both series do not include sporadic vi- 
brations because the duration was too short to accomplish 


measurement. 
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TABLE VI 
IMPRESSED HEATER VOLTAGE 
(Body of table is volts) 


Run Series 
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TABLE VII 
Run Series 

D 


ELECTRIC HEATING CURRENT 
C 


(Body of table is amperes) 
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INITIAL AIR STREAM VELOCITY DATA 


(Body of table is velocity, feet per second) 


Run Series 
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TABLE IX 
Run Series 
D 


(Body of table is degrees Fahrenheit) 
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INITIAL AIR STREAM TEMPERATURE DATA 
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Run 
No. 
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TABL= X 


SiRIES A RUNS 


Xy 


in. 


O ADNAN OO UII SI ON PO NINO 


J ENOL NS O69 F OR OBNO NIN RPO 


THELPRRATURE DISTRISUTION DéATa, 


(Body of table is ty-to, F) 


Run No. 

2 3 4 5 6 7 S 
Vea 5 (SRO (13.6 lay! 1ea2- eee Via 
14.8 17.3 15.4 15.5 15.3 13.8 15.8 
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17) VeOul (1s e WSet lyao sro ieee 
I7et 20%6 We.7 Vas Lea Wees ar 
Were GeO. 7  1s.8° Vous oe Dees ace 
18.4 ae att lii.o 12.4 16.5 We 
iesO SS 20.+ le. o) Wed eee a 13.2 
7.9 “20:2 Loss -11e.2 "ewe ee ies o 
17.3. 19.26 D729 “17.7 179 W6s1 eG 
$Oeo) Ven eel al ie eee ee 
lo,5 Werl “lols <ieve.- 16e7% “ies ode. 
15.0 17.1 15.5 15.5 15.9 14.5 16.2 
ae 7 25.7 Visa ae 14.8 17.6 15.2 
12.1 13.4 12.3 12.4 12.8 12.1 13.5 

G0 “10.2 “9.00 Ge yor? “676 Ste.e 

Run No. 

ipl 12 ik 14 15 16 
16: el ee eee ee 
16.9" 15.59 13.2 Teese eng ie, 

13.0 7 1 le tee 17 
16-5 17.4 Vigee low ee; 15.9 
16. ¢ “1¢.5- aleey Toe ea! en / 
19.2 — yy prolate SN) 
Voee (eet ade ee sites “eee 
19-5 (2.9 (17. «Tose ee ae 
19.9 19.3 18.1 165.4 185.4 ie 
19.7 “19,4. de levteateee “ee 

19.7 19.3 17.8 16.4 16.2 16.3 
19.5 19,0) Sy. ee ee Oe 
1.2 (1S 7 i oe ee ee ye 
16.7 “18.3 G70 ee awe 7c 
17.0 17.0 a6st ioe lon 
13.3 12. Selene Bemmeeewe le 


NM OCO™N PWN FON ONOUITIWIU? 10 08 





3) 


TABLE X, Continued 


THMPHZRATURE DISTRIBUTION DATA, SERIES B RUNS 


(Body of table is ty-to» F) 


Xy 


i. 


Run No. 

iL 2 3 5 6 7 & 
eu,4 22,7 21.7 18.9 17.6 16.2 15.8 14.5 
26.6 e4.6 23.8 20.5 19.4 17.6 17.4 16.2 
Pees 26.6 25.4 21.6 20.9 19,2 19.1 27.7 
Bem 27e4 “26,e 23,2 2le7 20n1 42000, Wee 
Boes> 27,6 26.7 23.7 22.5 20.56 20.6 19%; 
21.1 26.0 27.4 24.1 22.9 21.2 20.9 19.8 
31.2 28.0 er 24.1 22.9 21.4 21.1 19.9 
21.1 28.0 ag eee 2 sel eit ions eee 
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26.9 2 ¥ re C289 22.1 COnf “2009 Veo 
Be e>.4 24.4 (21.8 21.2 1969 49.9 Vige7 
Bo,0 23.9 23.2 ar 20.4 16.9 19.6 19.0 
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Wee, 15.3 17s Oso 2Se9 15-2 15.7 “525 

= 12.9 2.5 12.07 Viet (es “il foe 
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10 11 12 13 14 15 16 
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TASLE A, Continued 
TEMPERATURE DISTRISUTION DATA, SERIES C RUNS 


(Body of table is tet yt) 


Xy Run No. 
oi. #861 2 3 4 5 6 1 g 
O 25.4 ea4u,4 22,1 21.5 22.0 21.1 21.2 20.3 
Oe 27.2 25.9 23.6 ot 23.6 22.f 23.1 22.3 
RH 29.0 27.7 25.3 24.6 25.6 24.4 25,0 24.1 
oe 30.0 26.9 26.2 a 260.8 25.6 26.2 25,4 
S “eat 20 eee], ee. 27.0 20.4 27.1 26.1 
fo 671.4 30.3 27.5 27.2 28.6 27.4 28.4 27.5 
mM. 70.9 30.0 27.0 26.8 28.5 27.2 26.4 27.7 
me 30,4 29,4 26.4 26.3 24.2 2 19 26.1 21. f 
Mm 29.2 28.7% 25.5 25.7 27.4 26.4 27.8 27. 
me e7.4 27.1 24.2 24.% 26.1 25.3% 26.8 26.6 
2c0 24.8 24.5 22.0 22.2 23.9 23.2 24,8 25.0 
Pe 0.4 20.0 17.9 18.2 19.6 19.0 20.5 21.0 
x» Run No. 
in. 10 skal 12 13 14 15 
O 20.7 19.8 19.4 19.5 16.1 17.9 
Mmmeee.6 621.5 «621.0 «621.0 19.2 19.2 
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TABLE X, Continued 
TEMPERATURE DISTRIBUTION DATA, SERIES D RUNS 


(Body of table is tanto» F) 


Run No. 
2 3 4 5 6 7 & 
O- Gese9) Cbs Veil 2 ©2220) =2inG eeonen 27 
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TABLE X, Continued 
TEMPERATURE DISTRIBUTION DATA, SERIES E RUNS 


(Body of table is t,-t,, F) 


Run No. 

2 3 yy 5 6 ii & 
BO 2155) 1925 1748 <6e6. BEG Goa aw 
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CHAPTER VII 
PHYSICAL PROPERTIES AND OTHER CONSTANTS 

Physical properties and other constants are required 
in order to rencer the observed data into a more useful 
form. It was noted previously that the relative humidity 
of the air was a maximum of thirty-five per cent at the 
blower inlet. Such a relative humidity is small enough so 
that the air may be considered to be dry for present pur- 
poses. As a consequence, the properties of dry air were 
used and were obtained from the NBS-NACA Tables of Thermal 
Properties of Gases. These properties include iene enle 
viscosity (26), thermal conductivity (27), and Prandtl 
number (28). 

An estimate of the thermal conductivity of the mater- 
jals of construction is required also. The value chosen 
for the thermal conductivity of both inner and outer brass 
tubes is 59 B hen ae according to Jakob (5). The 
value chosen for the thermal conductivity of the heat- 
resisting steel outer sheath of the Calrod heater is 
1 pvp 


15 Bhr f , according to Schack (29). The value cho- 


sen for the magnesium oxide insulation in the Calrod heater 


is 0.2 B nr ft tpt 


» assumed to be equivalent to mica, 
according to Jakob (5). 

Some of the heat transferred from the surface of the 
heat transfer element occurs by radiation to the surrounc- 


ings. As the surface was bright chronium plated and then 


brightly polished, the value chosen for the hemispherical 
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t energy conversion factor ised was 3.413 Bae 
att-hour, according to Keenan and Keyes (25). 
The computed areas of the various parts are: 
Part Area of cross-section, £t°. 


iter tube, excluding holes 1.238 x 1077 


nner tube, 4 
excluding thinned section 1.914 x 10 


Salrod outer sheath and coil 2.09 x 1077 


lrod insulation ano es 100% 
The computed outer circumference of the outer tube is 


0.16 34 1G. 











CHAPTER VIII 
METHOD OF CALCULATION 
The method of calculating the mean surface coefficient 
of heat transfer is based upon a heet balance and a defin- 
ing equation. A heat balance may be written as: 
a a Ge a te dowon 
Q = Gy- Wn - Weg - Gy - VIII-1 
Substitution from the defining equation, Equation II-%3, 


yields: 


=o) Sue ; Vill=e 


in which each term on the right must be evaluated. 
The term qs is the electrical heat energy released in 


tne thermal length x If the heating coil length is My 


i’ 
then the heat release per unit of length is (q,-/1) and the 
heat release for a length x; is (x ,q),/M) for uniform heat 
release. <A correction factor may be applied if the heat 
release is non-uniform. With such a factor, the heat re- 
mease becomes: 
Qy = srccc--- : VIII-2 
M 

Values of q.. are given in Table XI. The data in the 

table come from multiplying the product of corresponding 


values in Tables VII and VIII, respectively, by the conver- 


sion factor 3.413 B per watt—hour. 


6g 
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A coordinate system is chosen to describe other terms, 
The noseoiece joint is selected as the origin and the vari- 
able x is chosen to increase in a direction downstream from 
the joint along the heat transfer element. 

The term 96, is defined by Equation II-5 and may be 
written as: 

pl 
@,, = (1/x,) JO dx. VIII-4 
0 

The term a is the heat lost to the nosepiece and may 

Perwritten as: 
de 
Gee a KYA e a . VIII-5 
The term gg is the heat passing along the heat transfer 


element downstream from the section located by x; and may 


be written as: 


do 
Oe > Re ; VIII-6 
dx X = xX, 


The term q, accounts for the heat lost to the surround- 
ings by radiation from the surface and may be written as: 


q h,Cx,9,, ; Viti" 


r 
Data for 9 are not yet available. However, the data 


from Table X may serve for estimation purposes. Heat con- 


duction through a round cylindrical wall may be expressed as: 


Gg = «= eT best <aee--- as 
Lie ee) 


a ee 





gel 


ier 0 OS ae eee VIlles 


Fema fixed construction, Equation VIII-S shows thee 
the temoerature difference is a maximum when qo 1s a max- 
imum. The maximum heat input in any run was 199.5 B hrvt, 
found in Table XI. The ratio (ro/r,) may be written as 
(0.624/0.384). The value for k, may de 59 8 hr tt tr}, 
and the value for L, may be (gee oy eo epee 

Solution of Equation VIII-& may be accomplished with 
the above values if it is assumed that all the heat gener- 
ated passes through the wall. The solution yields t,-to 
to be 0.099 F as a maximum. On the other hand, the place 
of measurim: the temperature was at half the wall thickness 
So that the temperature difference between the surface and 
the place of measurement is 0.05 F, approximately. 

There is close geometrical similarity between the pres- 
ent instance and that of Jakob and Dow (2), and their study 
Oi the temperature field near the holes in the outer tube 
wall may be applicable. Their ratio of outside diameter to 
inside diameter is (1.3/0.95) equal to 1.37, the present 
one is (0.624/0. 354) equal to 1.62; their ratio of wall 
thickness to hole height is (0.175/0.08) equal to 2.19, the 
oresent one is (0.120/0.0625) equal to 1.92. 

It is considered that aerodynamic heating may not in- 
fluence the temperature measurements. The relative velocity 


between the surface and the thermocouples used for measuring 


air stream temperature is zero. Both items have the same 
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general shape, that of a round cylindrical probe, Although 
the heating effect may have been as great as 1 F in some 
runs, it is considered that the effect on each item is the 
same and that the effect applied to the difference between 
the measurements is negligible. 

AS a consequence of the discussions above, the data of 
Table X are assumed to be applicable as measurements for 
Surface temperature, The surface temperature may be found 
by adding corresponding values in Tables IX and X, respect- 
ively. 

The values in Table X were plotted as a function of 
thermal length for all runs as illustrated in part in Fig- 
ures S12, inclusive, to fair or smooth the data, 

The irregularities in the data of Table X may be due 
to (a) hot spots, (bv) air stream temperature stratification, 
(c) variations in thermocouple placement, and (da) the data 
are the difference between two magnitudes. 

Most irregularities appear at high energy input rates, 
It is considered that hot spots are the probable cause of 
irregularities, because of the method of construction. The 
parts may touch at some places and may be entirely separa- 
ted at others. 

Stratification in the air stream was found to occur, 
the extent depending uoon atmospheric conditions. This 
effect was the underlying reason for establishing tempera- 
ture variation limits as recorded on the Speedomax, dis- 


cussed in Chapter VI. 





ec 


Variation in thermocouple location may have been one- 
eighth of an inch for it is assumed that placement was 
mace to plus or minus one-sixteenth of an inch. An overall 
variation of placement of one-eighth of an inch could cause 
an irregularity of O.1 F in the neighborhood of high temper- 
ature gradients along the element. 

If it is assumed that the primary temperature determin- 
ations are within plus or minus 0.1 F then a difference of 
two such measurements could yield an irregularity of 0.2 F. 

The smoothed data are found in Table XIII. The data 
from this table are used in subsequent calculations. Sim- 
ilarly, smoothed data for surface temperature may be ob- 
tained by adding corresponding values in Tables IX and XIII, 
respectively. 

The mean temperature difference is found from Equation 
VIII-/7 by applying Simpson's Rule for integration. Use of 
the rule affords vorogressive integration with ease and per- 
mits choosing thermal lengths in sequence at intervals of 
four incnes. Several spot checks were made on the results 
of integration by using the data of Tables X and XIII, res- 
pectively. The difference between results did not exceed 


one per cent for the trial checks that were made. 


de ad 
The derivatives |{-- and { -- are approximated 
ax / x=0 ax X=X, 


by finite differences. The former derivative is approxima- 


; -l 
ted by (05-0,)/(X5-X) equal to 6(0,-8,) F ft -, hence 


il 


d, = 6k,A(@5-85) B hr VIII-9 
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The latter derivative is aporoximated by assuming that 
the best value is the average slove of the two secants 
joining the points (x,_9, 94_5), (xy, 9,4) and (x,, 94), 
(X54) Os4o)s respectively. The slope of the one secant is 
(0,-9,_5)/(X4-X;_9) equal to 6(0,-0,_5) F ft}. Similarly, 
the slope of the other secant is 6(03; 40-9; ) F ft, The 
average of the slopes is 3(0;40-9;_59) F ft 7, hence 

de = —3k,A(@;,0-0; 5) Bahr? . VIII-10 

The evaluation of the product (k,A) is based upon the 
assumption that the temperature gradients along the outer 
tube, inner tube, and heater are the same although the tem- 
perature level of each may be different. It may be supposed, 
as an example, that the parts are centered. The air space 
between the inner and outer tubes is 0.004 in. and is 0.005 
in. between the inner tube and heater. The thermal conduct- 


-l,,-lp-l 


ivity of the air may be 0.0156 B hr ft “F -at 100°F and 


the maximum and winimum heat release rates from Table XI are 
mo.5 and 23.6 B nrvt, respectively. With the aid of Equa- 


tion VIII-&S, tnere results; 


Power Input tj- to, F 
Beeer hr O.004 in. space 0.005 in. space 
199+) 1520 Ze. 
eoy ie eno 


The above data are indicative of the respective temperature 


levels. 
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The calculation of the product (k,A) is as: 


Part k A kA Per cens 

of 

Bhrvgt tr) ft? B ft hrtr-l ss Total 

Outer tube 59 1.230) % NOer 73.1 x 10-2 63.6 
Inner tube 59 Toma yoy, wie lon 12.9 
Calrod sheath 15 2.09 x 107 3.1 x 1072 3.5 
Calrod insulation 0.2 3.34 x 107 0.06x 1077 a 
Total $7.5 x 1072 100.0 


This result shows that 83.6 per cent of the heat con- 
ducted along the heat transfer element is conducted by the 
outer tube, if the temperature gradients are the same for 
each part. 

An estimate of the consequences of the assumption is 
of interest. At the nosepiece joint, for example, the tem- 
perature gradient in the outer tube was limited to 1 F per 
inch or le F per foot, approximately. It may be Supposed, 
for example, that the gradients along the inner tube and 
the heater are in the proportion 1:2:4. For these condit-~ 
ions, the ratio of the new kA product to that computed 
above is 1.c4%. With a heat input rate of 199.5 B nro, the 
heat input to a length of four inches is 25.e B hr-+, The 
loss through the nosepiece in the one case is (12 x 0.0875) 
equal to 1.05 B nro, or 4.17 per cent of the heat input. 
In the other case, the loss is (1.05 x 1.24) equal to 1.30 


B nr-t 


, or 5.15 per cent of the heat input. Thus the loss 
error is 0.95 per cent of the heat input. 
The same calculation with the minimum heat input, 23.6 


B hr +, yields 35.3 per cent and 43.6 per cent of the heat 
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input, respectively, and the loss error becomes &.3 per 
een. 

The loss error is seen to decrease as the length is in- 
creased. At those places where the temperature gradient 
along the outer tube approaches zero, the basic assumption 
becomes more valid. As a consequence of the above discuss- 
ions, the temperature gradients along the members are assumed 
to be the same as that along the outer tube. 

The correction factor Iie may be estimated. Measurement 
of the heat input to a heating element is seen to yield the 
total value and no information as to release distribution 
can be obtained from such a measurement. Thus two coils 
having the same electric resistance may be wound with uni- 
form wires so that the wires are bunched together in one 
coil and evenly spaced on the other; or, a uniform wire may 
be used on one coil and a tapered wire on the other. The 
measured heat input rate would be the same for all. 

Let it be assumed that the temperature distribution 
along a coil of length Mis given as 9 = a, - box. The heat 
ieamsierred is; 

Mi 


se hoe dx , and for constant h, 


and this is also the measured heat input, that is, q = q,. 


lied DG is a small number, h may be assumed to be constant so 





= OO, DREN ee Per Pre Cig eee Aicow = ) «fe. (esas. 
Siubhetitutian for aq vielde: 
Substitution for q yields: 
ions = erties Peet arr ee ge ee 3 sa Ren a 


2 
boxy 
SANG Eee eee 


ya 
The assumed heat input from input measurements may be q, 
Cx 58 


m? Which is less than q;, becoming equal to it at 


The factor F, may be defined by the ratio (q,/q,). 





BS ee VIII-ll 


It is recalled that the Calrod heater tests could be | 


represented by the equations: 
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Pn 


On = 160.4 ~ 0.19 2. 


367.1 - 0.43 z, and 


Figure / shows that the equations may be extended to the 
left beyond the middle twenty-four inches (z = 12 to 36 in.), 
the basis originally chosen. 

The relation between 2 and X, is given by X=) 2 =a ee 


By substitution, there results: 


oe 362.8 - 0.43 X, ir 


On = 160.4 + 0.19 X, = a, - bpx. 


The ratio (o./2a,) may be computed from these, for (b,/2a,) 


a,- Ox, and 


is (0.43/362.8) equal to 0.000593 and (b,,/2a.,) is 
(0.19/160.40 equal to 0.000599, with 0.000596 as the aver- 
age, 
Substitution of this result into Equation VIII-11l 
yields: 
1 ~ 0.000596 x; 


where M is 31.75 inches and x, is in inches. 
The correction factors for the calculation of the mean 
surface coefficient for the thermal lengths used are: 
7, 4 g 12 16 20 el 2g 
F.1.017. 1.014 1.012 1.009 1.007 1.004 1.002 
The heat transferred by radiation from a small body en- 


closed by a larger one is expressed by Stefan-Bolzmann's law 


and by Kirchoff's law as: 
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yoy 

T r 

O17 Pes (:.) = ey \ B nr? 
100 100 


Gr = 
= h S(T - T,), 
from which: 
T VE /t 
0.174 € a 7 — 
100 100 
aa 8 nr fe ope. 
T.- T, 
VIII-12 


The order of magnitude of h,. 16 Of Inverest. ~ it may 
be supposed, for example, that T, = 572°R, T, = 540°R, and 
fat © = 0.07. With the aid of Equation VIII-12, h.. 1s 


1 wend 


found to be 0.088 B hr ft If it is supposed that 


T, = 536°R, T, = 5e5°R, with other values as before, the 
computation yields h,. to be 0.075 B nev tet7opl, These re- 
sults reoresent extreme values encountered in the experi- 
ments. 

Also, the experiments were conducted during the winter 
season so that the temperature of the surrounding room 
walls was less than that of the air stream. As a conse- 
quence of the above discussions, a constant value of hy. is 
used in all calculations and is chosen to be 0.1 B hr7+rtr7}, 

Equation VIJI-2 may now be written in the form of a 
Merking equation as: 

Fyq.(x5/31.75) - 0.525(@5- 05) + 0.263(0,4,5- 94_5) 
0.1634 (x,/12)0 
By Simpson's Rule, 9, = (2/3x;)(@o+ 400 + 2Oyt+°°* + 6,) for 


an integration interval of two inches. Substitution for 0, 








&0 


yields the final working equation as: 


a a Sa ne a ae a a a a AD ED, SD SS SUS SN | SS EE oe [aa -O.1, 
(Oo+ 4Oo+ COyte++ + O;) 
VIII-13 


eal 


let | ae Xy is in inches, and corres- 


wnere h, is in B hr 
ponding 9, are from Table XIII. 

A local surface coefficient of heat transfer may be 
computed also, The computation is based upon an approxima-— 
tion to Equation II-l by finite differences over a thermal 
length of two inches. 

Let three points as (Xs _o» 8.5)» (x55 Q.), and 
(X54) O40) be chosen. The heat gained at a section at 
X,-1 inches may be approximated by: 
= Ee ee) 2) 


-1 
~6k,A(0,-0, 5), Bhr™ . VITI-14 


Gc,0 


The heat lost at a section at x,+1 inches may be app- 


roximated by: 


Gde.2 > -6kA(0,.5-9,); B hr? VIII-15 
>) 
The electrical heat energy released may be given by: 
2 
q, = F,q,(e/M), B hr : VIII-16 


where Mis in inches. 
The heat transferred through the surface may be written 
as: 
q = (h,th)c@,(2/12), Bhr. ange 
The heat balance for the portion of element may be 


written as: 


Ie,o0f Ai = ee Q VIII-1& 
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Substitution into Equation VIII-18 from Equations 
WITI-14 to VIII-17, inclusive, yields: 
' 
: Fq)(12/M) + 36kA(0,,5-0, 5-20, ) 
Substitution of previously determined constants gives 
the final working equation as: 
t 
2. 305F,.d),-19. 28( 205-9, 45-03_9) _ nr-tpe-ep-t, 
at VIII-19 
1 
in which F, is yet to be determined. 


t 
The factor ee is defined and determined as: 


Ky + | 
hce dx 
x,-l 
i = Le San a en ea 
ehce 
2a 
On 
and , 
a -b xX 
ia eee a eee oe 


O w 
me Lee 
0 

l= = 
cay 


Substitution of the values previously determined yields: 


. 1-0.001192 x, 
Py = wee eean-----= VIII-20 
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The correction factors for determination of the local 
coefficient of heat transfer for the thermal lengths used 


are: 


Fy Seo ( 0.995" 02995 0.991 “Ol 9ee. von 9e5 0.95% 

It is recalled from Chapter III that the Reynolds 
humber offers a description of the flow configuration. The 
characteristic length is best chosen as the length upon 
which the boundary layer is formed whether heat is trans- 
ferred or not. Therefore, the characteristic length is cho- 
sen as the total length and the kinematic viscosity is cho- 
sen at the undisturbed air stream temperature, vo: 

The choice of the total length may be arbitrary in 
this instance. If the intersection of the longitudinal axis 
of the heat transfer element with the tangent plane at the 
nosetip is called a pole, then the starting length is de- 
fined as the distance along a meridian from the pole to the 
nosepiece joint. The thermal length is chosen as the dis- 
tance along the heat transfer element measured downstream 
from the nosepiece joint. With these definitions, the total 
length becomes the sum of the starting and thermal lengths. 

The starting lengths may be computed as: 

ABE « 


ee te, ee 
2 a 
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oo) Oia 
s=L + r( ---- = 1) 
i 2 
= Jo + 0.18, inches. VIII-~21 


The starting lengths used are given in Table XII. 
TABLE XII 


HYDRODYNAMIC STARTING LENGTHS 


Run Series Nosepiece length Starting length 

inches inches feet 
A ORE 0.49 0.0409 
B Dine ee 0.452 
C eco re 1.020 
D 18.06 15.24 1.520 
E e+.00 e4.1 ream OE Hy 
X Sees 5.43 0.452 
V sri, 5.43 0.452 


The concept of the Reynolds number as a fluid flow 
number may be applied to the Nusselt number. The Nusselt 
number becomes a heat transfer number. If the concept is 
continued, the total length in the Reynolds number is re- 
olaced by the thermal length in the Nusselt number. It may 
be noted that the characteristic length in both numbers be- 
comes the same when the starting length is zero. Therefore, 
the characteristic length in the Nusselt number is chosen 
as the thermal length. 

The choice of the physical properties of the fluid for 
the Nusselt number and the Prandtl number is based upon a 
film temperature. The film temperatures are defined as 
follows: 


For computations with h,, 


te = tot —-—-, F.’ VIII-22 
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For computations with h, 


a) 
t= tS 
f O 5 


ow VIII-23 


The procedure described shows how Nusselt number and 
Reynolds number may be computed, The general pattern of 
further calculations is to develop an equation in the form 
of Equation II-8&. 

The Prandtl number for air in the range of the exper- 
iments may be assumed to be constant. The extreme values 
are 0.710 at 65°F and 0.708 at S4°F, with an average of 
0.709. It is seen that the extreme variation of Prandtl 
number is in the order of one-third of one per cent. 

With this result, Equation II-S may be written as: 

ce Cees. VIII-24 

The procedure for subsequent calculation may be: 

(a) Establish regions of laminar and turbulent boundary 
layers, 

(b) Find a suitable exponent for Reynolds number for each 
Surface configuration, 

(c) Find a function go and 

(d) Establish C, and compare with values for other curva-— 


tures. 
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TABLE ZIII, Continued 


CORRECTED TEMPERATURE DISTRIBUTION DATA, SERIES B RUNS 


(Body of table is t.-t., F) 


Run No. 
ih e b 5 6 i & 9 
e444 22.7 21.7 18.9 17% 16.2 15258 14.5 Ie 
Gan> 25.3 235.9 20.5 19.6 17.9 17% rie 14 
meee 20.5 25,2 e.6 20.9 29.2 1 17.6 7 
29.5 27.3 26.1 22.9 21.8 20.1 20.0 18.6 1 
pees 27.6 26./ @23.{ @e.5 20e/ 20.6 Wels sag 
et 2f.l 27.2 24.1 22.9 21.1 21.2 19,7 Te 
peee 26.1 27.3 @F.e 23.) 27.8 2153 20.0 “ae 
pie 27.9 @7.3 @H.2 23,3 21.8 2155 2O,e “is 
30.5 27.6 26.9 23.9 23.0 21.3 21.5 20.3 18 
foo 27.1 26.3 23.5 e@2.6 21.1 Gl.3 etec “Le 
26.9 26.4 aan 22.9 22.1 20.6 20.9 20.0 18 
27.5 25.4 24.4 21.9 21.3 19.9 te 19.6 18 
25.7 23.9 23.0 20.6 20.1 ei 19.4 18.9 17 
egee 21.6 20.6 16.6 18:59 17:4 1850 17.6. 917 
mm 16)/ 7,0 15.6 Bees 155) 15.6 tone ets 
- N29 12.5 Tis Wiel 9e.> 1251 tons wie 
Run No. 

10 11 12 13 14 15 16 
fet 11.1 10.5 10.9 12.0 3259 1286 
12.9 13.1 12.5 12.9 13.9 14.9 15.1 
4,2 14.6 13.9 14.3 15.3 16.3 16.6 

15.2 15.7 15.0 15.3 16.3 17.4 17.6 
MeO Gr F559 sib o wy 0 Mtedy Peee 
mo.6 17.0 16.5 be.4 “27.6 “2657 18:7 
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15.0 15.3 14.9 15.0 16.5 17.1 17.6 

HO.7 10.5 11.0 10.9 24.9 13.2 14.1 
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TABLE XIII, Continued 
CORRECTED TEMPERATURE DISTRIBUTION DATA, SERIES C RUNS 


(Body of table is toto. F) 








Run No. 
1 2 3 2; 6 7 g 9 
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22 0621.4 20.0 18.2 18.2 19.86 19.0 20.5 21.0 22 
X, Run No. 

mam, 10 11 12 13 14 15 

O 20.7 19.86 19.4 19.5 18.1 17.9 
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6 26.9 25.6 24.8 25.1 23.1 22.6 
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o e o w e e@ 9 ® a e s ® 
OPN FOr MRE UIN™N 





e ® 


BUSBYBBRAS - 
OWN ONO LO Orn ONO © 


10 


DO 
aw 
1 ONO F OAXO OF OE 


O2RNO KO LO LO 010 CANT UI NM 
e e e e e © e e e @ 
Oa\n OVOR8 OR ON NO KO 


MO MP 1M 1M 11 PN) PO 


&8& 


TABLE XIII, Continued 


Mere ee 


NM MOM MO 1M fo 
NANDWAN FPWR 


Run No. 


22.0 
24.0 
2 

26.2 
2730 
24° 
ob. 8 
26.5 
26.0 


aie ae 


MN 
ae) 
MPrROMmFOW YW HA 


CORRECTED TEMPERATURE DISTRIBUTION DATA, SERIES 
(Body of table is t-to» F) 


ae) 
NM 
MMR AFRO 


D RUNS 





9 


TABLE XIII, Continued 
CORRECTED TEMPERATURE DISTRIBUTION DATA, SERIES E RUNS 


(Body of table is to-to> F) 
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Mm 29.1 26.0 24.1 21.7 19.6 17.9 15.0 17.4 18 
14 26.9 25.7 23.7 21.5 19.2 17.7 14.8 17-5 18 
i> 28.1 ee a1 Gls We7 75 BS 2 a 18 
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X, Run No. 
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4 19.1 19.3 18.0 a 
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TABLE XIII, Continued 





CORRECTED TEMPERATURE DISTRIBUTION DATA, SERIES X RUNS 


(Body of table is ty-to» F) 


X, Run No. 
in 1 2 3 4 5 6 7 gS 9 
O 18.86 19.6 18.3 17.0 15-9 15.2 14.9 11.4 11.5 
2 20.7 2ie5, @0.5 16.8 17.6 17,2 16.9 13 1 2 
wh 8622.1 22.9 21.9 20.1 16.9 185 18.4 14.3 14, 
6 23.1 23.9 22.8 21.1 20.0 19.5 19.4 15.2 15.5 
6 23.9 24.6 23.4 21.7 20.7 20.2 20.1 15.9 16.2 
ie 29.5 25.1 23.7 22.2 2162 20.7 20.6 16.3 16.7 
12 24,9 25.4 23.9 22.4 21.4 20.9 20.9 16.7 16.9 
me e540 25.8 28.@ 22.4 21.5 21.0 21.1 &6.6 17.0 
16 24,9 25.2 23.86 22.3 21.4 21.0 21.0 16.9 17.0 
16 24.5 2h.s 23.5 21.9 21.2 20.7 20.7 16.8 17.0 
20 23.9 24.2 22.9 21.3 20.7 20.2 20.4 16.7 16.9 
Mm 23.1 23.4 2212 20.5 20.1 1946 19.8 16.5 16. 
eh 21.9 22.2 21.1 19 29.1, 18.7 19.0 16.1 16.3 
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TABLE XIII, Continued 
CORRECTED TEMPERATURE DISTRIBUTION DATa, SERIES V RUNS 
(Body of table is tite» F) 


Xn Run No. 
in 1 2 3 5 6 7 & 
O° 26.1 22.0 20.6 18.1 14.5 13.4 13.1 23.1 
@ 28.3 23.7 22.7 19.9 16.7 15.8 15466 15.8 
4 30.2 =o 24.3 21.5 M.e Fe 17.1 abe 
6 31.6 26.0 25.5 22.9 19.8 18.4 18.3 18, 
6 32.6 26.5 26.5 24.0 20.6 19.2 19.1 19.4 
10 33.2 26.8 27.0 24.9 21.5 19.8 19.7 19.9 
M2.) 3302 20.6 27.2 25m 22.1 20.28 20a 20,2 
14 229 26.6 — 25.6 22.3 20.5 20.3 20.4 
16 32.4 26.2 26.9 25.9 22.5 20.7 20. rg 
16 31.7 25.6 26. 233 22.6 20.9 20.5 20. 
20 20. f 25.0 26.0 26.0 22.8 21.0 20.6 20.6 
220 289.4 24.1 25.2 ere 22.9 21.0 20.5 ae 
2h 27.6 22.8 24.0 24.6 22.7 20.9 20.3% 20. 
26 625.1 «9220.8 —¢ 26.0 21.6 20,4 19.7 20.0 
26 21.6 18.0 19.4 20.3 19.7 18.7 18.2 18.7 
30 «416.0 13.4 14.5 15.1 15.3 14.5 14.4 15.0 








CHAPTER IX 
EXPERIMENTAL RESULTS 

The preliminary experimental results are developed by 
expressing Nusselt number in terms of h,. It may be noted 
that Equations VIII-13 and VIII-19 offer a choice between 
h and h, for calculations. 

It was shown in Chapter III that all theoretical re- 
sults and one experimental result for determining the func- 
tion g™ in which h is to be used have been presented. One 
experimental result is available in which h, is used. This 
would indicate that calculations with h may be preferable, 
for better comparison may be possible. 

Closer examination of Equation VIII-19 shows that 
errors of approximation may be appreciable. The equivalent 
of MeO / ax") is required and the attainment of an accurate 
approximation is a formidable obstacle. The factor F. is 
of importance and should be known accurately. These consid- 
erations show that the calculation of h may be of acceptable 
status only near the place where the temperature gradient 
mong the outer tube is near zero. 

Examination of Equation VIII-13 shows that the equiv- 
alent of (dQ@/dx) is required as well as the factor Fl. How- 
ever, it has been shown how the effects of errors of approx- 
imation may decrease as thermal length is increased. Asa 
consequence, h, is chosen for calculations, 


The experimental results are shown graphically in 
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Figures 13 to e3, inclusive. Nusselt number is plotted as 

a function of Reynolds number with logarithmic coordinates, 
The results illustrated in Figures 13 to 20, inclusive, 

were obtained with the air stream flowing parallel to the 
longitudinal axis of the heat transfer element. It is noted 
that laminar boundary layers may have been obtained for 


(Npe)y up to 10° 


and that turbulent boundary layers may have 
been obtained for (N,.), at 3.x 109. It is also noted that 
the exponent for Reynolds number is verified to be 0.5 in 
the laminar range and to be 0.8 in the turbulent range; the 
lines drawn in the Figures have corresponding slopes. 

Each of the Figures is drawn for a specific thermal 
length. Comparing Figures 13 to 19, inclusive, it is noted 
that at a given Reynolds number the values for the different 
starting lengths (Series A, B, C, D, and E,) come closer to 
each other as the ratio s/L becomes smaller. This effect is 
noted for both laminar and turbulent flows. 

The effect of oblique cross flow, with an angle of three 
degrees between the stream and element axes, is partially 
shown in Figure 21. The starting and thermal lengths chosen 
for the illustration are 5.43 in. and 16 in., respectively. 
At Reynolds numbers greater than 5 x 107, the slope of a line 
drawn through the data corresponds to an exponent of unity. 
At Reynolds numbers less than ie & Oey the slope corresponds 
to an exponent of 0.5/. 

The behavior at high Reynolds number may have been an- 


ticipated. The surface coefficient of heat transfer in 
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rectangular cross flow is independent of the length for 
given fluid velocity and properties, excluding end effects. 
If the distance along the surface is chosen as character-— 
istic length for both Nusselt number and Reynolds number 
then the variation of Nusselt number with Reynolds number 
should become linear; that is, the exponent for Reynolds 
number should be unity. 

A possible conclusion, therefore, is that data in the 
turbulent range suggesting an exponent greater than 0.8% 
for Reynolds number with respect to Nusselt number may have 
been influenced by oblique cross flow, 

The behavior at lower Reynolds numbers, approximating 
the transition region, is unexplained and is of no import- 
ance in the present work. 

Comparison of Figures 13 to 20, inclusive, with Figure 
Cl shows that centering of the heat transfer element in the 
air stream was accomplished with reasonable effectiveness, 
The exponent of Reynolds number remains at 0.8 in the runs 
with parallel flow. 

The effect of vibrations of the heat transfer element 
is partially illustrated in Figure ee. It is seen that the 
exponent for Reynolds number remains at 0.& for the vibra- 
tions measured in the turbulent range. The effect of vi- 
brations in the laminar range is not clear. 

Composite results comparing vibration data, cross 
flow, and parallel flow are illustrated in Figures 23 to 


e&, inclusive. At high Reynolds numbers, the exponent for 
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cross flow is seen to remain at unity and the exponent for 
vibration data remains at 0.8. It is noted further that 
the effect of vibration may not be appreciable. 

The effect of vibrations in the turbulent range may 
be examined quantitatively. The data of each run may be 
expressed as: 

(Niu)x = Fhe 
for the turbulent range. Rearrangement yields: 

C. = (Ng) / (pee? 1X-1 
in which C, may assume different values for different sur- 
face configurations. 

The computations for the Series V runs are illustrated 
in Table XIV. The data for the table were selected so that 
Reynolds number is equal to or greater than 5 x 10°. 

A similar computation for the Series B runs may be 
usec for comparison purposes for the surface configurations 
are the same for both series. The results are summarized 
in Table XV. 

TABLE XV 
COMPARISON OF VIBRATION AND PARALLEL FLOW RESULTS 


8, ie, Average Value of Cg Difference, per cent 
ia. in. Series V Series B of Series B Values 
5.43 8 0.018 0.0191 3.1 
5.43 12 0,021 0.0218 1.6 
5.43 16 0.0234 0.0238 oe 
5.43 20 0.0249 0.0253 i 
5.43 24 0.0261 0.0263 0.8 
5.43 28 0.0267 0.0269 0.7 


These results suggest that the effects of vibration 
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are within the limits of experimental and approximation 
error and that results with respect to the vibrations en-— 
countered may not be significant, 

Two conclusions are possible with respect to vibra- 
tions. It is possible that the combination of different 
frequencies and the respective amplitudes is such that the 
net result is the same for the present experiments. On the 
other hand, it is possible that the kinds of vibrations en- 
countered were of negligible effect in comparison with a 
system that does not vibrate, 

Eckert (10) has stated that the greatest resistance to 
heat transfer by forced convection with turbulent boundary 
layers is due to the laminar sublayer, and that this layer 
is very thin with respect to the overall boundary layer 
thickness. It may be concluded, therefore, that vibrations 
have little effect on the heat transfer unless they are of 
such a nature as to interfere with the laminar sublayer., 

Examination of the data coupled with the premise of 
Eckert leads to the conclusion that the effects of the vi- 
brations encountered in the turbulent range are negligible 
with respect to a steady system. In addition, the conclu- 
sion is reached that vibrations may have a greater effect 
in the laminar range because the whole boundary layer is 
postulated to be laminar. 

Illustration (C) of Figure 1 may serve for illustrat- 
ing these conclusions, If it is supposed that part of the 


boundary layer to the right of L, is removed or washed 
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away, yet allowing the laminar sublayer to remain intact, 
then appreciable changes in the heat transfer may not be 
expected. On the other hand, if a part of the boundary 
layer to the left of be 1s washed away, then appreciable 
changes in the heat transfer may be expected, 

An additional effect may occur when vibrating and 
steady systems are compared. In the case of flow in pipes, 
for examole, it is known that laminar flow may occur for 
large Reynolds numbers. If such a flow is disturbed, then 
the flow becomes turbulent and remains so. Further, if a 
flow is occuring for Reynolds number less than the lower 
critical, then disturbances may introduce a turbulent flow 
which returns to ani aae upon cessation of the disturbance. 
Thus the effect of vibration may be to change the lower 
critical Reynolds number, an amount depending upon the na- 
ture of the vibrations. 

The supporting results of calculations for the Figures 


presented in this chapter are given in Tables XVI to XXII, 


inclusive. 
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TABLE XIV 


CALCULATION OF VIBRATION DATA, SERIES V RUNS 


Run No. 
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18. 00 
10.06 


2670 
2.65 


2.61 


28 
5.96 
Me 7 
1100 
2,64 
&. 33 
5645 
1457 
2.07 

11.66 
7.18 
1985 
2.64 
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EXPERIMENTAL RESULTS FOR THERMAL LENGTH OF FOUR INCHES 


FIGURE 13. 
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EXPERIMENTAL RESULTS FOR THERMAL LENGTH OF EIGHT INCHES 


FIGURE 14, 
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FIGURE 20. EXPERIMENTAL RESULTS FOR STARTING LENGTH OF 
O.49 INCHES, LAMINAR FLOW 








107 








NUSSELT NUMBER, (Naw) x 





10? 2 3 & 5 678 10° 2 


REYNOLDS NUMBER, ( Ne eL 


FIGURE 21. EXPERIMENTAL RESULTS FOR OBLIQUE CROSS FLOW FOR 
THERMAL LENGTH OF SIXTEEN INCHES. SERIES X RUNS. 
ELEMENT AXIS AT 3° TO JET AXIS. 
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FIGURE 22, #XPERIMENTAL RESULTS FOR VIBRATION FOR THERMAL 
LENGTH OF SIXTEEN INCHES, SERIES V RUNS. 
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CHAPTER X 
CORRELATION FOR THE MEAN SURFACE COEFFICIENT 

The computation shown in Table XIV suggests a method 
to determine the effect of surface configuration. The sub- 
sequent calculations are limited to the data for turbulent 
boundary layers. The data available for laminar boundary 
layers are not sufficient for detailed conclusions and the 
effect of vibrations in the laminar region was not explored. 

Data for correlations are, therefore, only chosen from 
Tables XVI to XX, inclusive, so that Reynolds number is 
equal to or greater than 5 x Ho? The procedure for compu- 
tation was the same as that outlined in Table XIV. 

A summary of the results of computation is given in 
Table XXIII. These results suggest a relationship between 
C of Equation IX-l, and the dimensionless ratio of thermal 
length to total length, (x/L). The existence of such a re- 
lationship is shown in Figure 29. 

It is noted that the results shown in Figure 29 offer 
a strong correlation. The range of thermal lengths consid- 
ered is in the proportion /:1, the range of starting lengths 
is in the proportion 55:1, and the range of the ratio (x/L) 
is in the proportion /:l. 

An equation for the line shown in Figure 29 was ob- 
tained from the data of Table XXIII by the method of least 
squares. The result of the computation yields: 


Co = 0.0307 tama X-1 


125 





126 


By definition, Equation IX-] 
: — 0.8 
oon CN) Mp) y 
therefore: 


(Ny, ¢ = 0-0307 (x/1)°°72( m2 °F K-2 


fee 
This result may be put into the form of Equation II-§& 
as: 
(Mya) t oa (Ni) x (L/x) 
0.0307 C717 aie 6 for air. X-3 


i 


If it is assumed that the Prandtl number for the air con- 
ditions as encountered in these experiments has the value 
0.709 and an exponent of 1/3 for the Prandtl number, then 
0.0 0.8 al 
(Myy)q = 0-034 (x/)°°97 (p)2*8( 24? 


is a general equation. a 
By analogy with Equation II-8, it is seen that: 
Cc, = 0.034 
g™ = er, go Xa5 


The results of correlation are illustrated in Figures 
30 to 37, inclusive. For purposes of illustration, a dimen- 


Sionless parameter is defined as: 


OF 8) 
Y = (Ny) ,/(x/L)°°7*, X-6 
and for the turbulent region: 
x 0.8 


Figures 30 to 36, inclusive, may be compared with Fig- 


ures 13 to 19, inclusive, as follows: 
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Before correlation After correlation 
Figure a compares with Figure 30 
i 1 t T | iT) 31 
i" 5 " in " 32 
(1 1 ? it it " S 
i 1 7 1 i " 3 
HH 1 § " it W a 
it L 9 } | it 11 3 


Correlation of all data is illustrated in Figure 27. 
It may be noted that correlation is strong for Reynolds 
number greater than 5 x 109. The result of Jakob and Dow 
(2) for the laminar range is included for comparison pur- 
poses. Those points lying below the Jakob and Dow line are 
based on low velocities which cannot be measured exactly. 

It is of interest to compare the results for the func- 
tion %” with another observer. A direct comparison may be 
made with Equation III-1, that of Jakob and Dow (2). The 
present result in the form given in Equation X-5 is conven- 
ient and the comparison is made neglecting the small effect 
of the different definition of starting and total lengths 
as used in Equation III-1. The comparison is made in Table 
XXIV. 

Equation X-5 was obtained with (x/L) having extreme 
values of 0.141 and 0.983. Equation III-1 was obtained with 
(x/L) having extreme values of 0.394 and 0.899. The data 
given in Table XXIV for (x/L) in the interval 0.4 to 0.9 
show a maximum deviation in the comparison of 1.5 per cent. 
This is a most satisfactory agreement, considering that the 
influence of surface curvature (to be considered in Chapter 


XII) cannot change much in the longitudinal direction. 





128 


TABLE XXIV 


COMPARISON OF SURFACE CONFIGURATION FUNCTION VALUES 


(x/L) 


OO ON AN Fu MH O 


MrPOOOOOCO0O0C00 


(s/L) 


ODWDVOOO0O00F 
OF MU FUT ON OO © 


Equ. X-5 Equ. 


1.400 
1.299 
1.216 
1.150 
1.09% 
1.059 
1.032 
1.015 
1.005 
relenl 
1.000 


IiI-1 Dirt. per cent, 


based on Equ. X-5 


Oh KER ORPWUIVI 


OWOUIIWN ROH PH & 
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TABLE XXIII 
CALCULATED RESULTS FOR CORRELATION PURPOSES 


Data selected for Reynolds number equal to or greater than 


5 x 102 
Run Series 5 in. x, in. x Lp Average value of CO, 
A 0.49 & 0.942 0.0280 
12 0.961 0.0257 
16 0.970 0.0290 
20 0.976 0.0290 
ek 0.980 0.0295 
28 0.983 0.0305 
B 5.4% & OSS 0.0192 
12 0.689 0.0218 
16 0.746 0.0238 
20 0. 7&6 0.0253 
24 0.815 0.0263 
2s 0.837 0.0269 
C 12.24 uy 0.246 0.00835 
& ON te ae 
12 O. 2? 0.01 
16 o- 0.018 
20 0.620 0.0204 
D 18.24 Ly Oneal 0.00650 
S 0.305 mre 
12 0.39 . 
16 5 ee 0.0158 
E e4u,1% 4 0.141 0.00510 
Sg 0,29 0.00870 
12 0.332 0.0115 
16 O. 396 0.0137 
8 
(yy) x/ Re) 
= gs + x, inches 
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CHAPTER XI 
THE LOCAL SURFACE COEFFICIENT 
Previous discussions have been concerned with the mean 
surface coefficient of heat transfer. An expression for 
the local surface coefficient may be developed from the 
foregoing results. 


Equation II-4+ may be written as: 
x 
DP m* = {ne ax Xian 
O 
Differentiation of Equation XI-1 yields: 


ad a x 
-- (h,9,%) = -- [re dx 


ax ax O 
and, d 
ho = -- (h, Ox) 
ax 
OT, On hynx den XOm dh 
h = hywn- + eee eee +t ere one XI-2 
.) -) ax 6 ax 


Equation JI-5 may be written as: 
x 
Onx = fe dim ~ XI-3 
O ) 


Differentiation of Equation XI-3 yields: 


do 
@=0,+x--- 
ax 
or de @- 9 
an Seaeeule XI-4 


Equation X-e2 may be written as: 


j Xx ~ V (gs er x) _ 
(Nyudx = So (=% & —-o——— 


Ss +x 
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z 
Ck /v ee ee Ne = 
x V; Ss +x 


which is a function of x only. Terms may be combined as: 


= oe 5 =k (s + x) tae e+ XI-6 


n 
dh,, ' x . Vo(s + x) Ll dk (m-1)k (n-m)k 
--- = Cy{ ----- | | --------- ~ -- + ---9-- + -------- 
dx s+ x Vo x dx x x(e + x) 
¢ 
dh, Cok ( x )" e& + = yf. ak (m-1) (n-m) 
OF, mee = mem f eemen— ff ae - ae +t nn + -—-—— 
dx x \s+x VW 7 k ax x (s+x) 
XI-/ 
Substitution from Equation XI-5 into Equation 
XI-/ yields: 
ah, 1 dk (m=-1) (n-m) 
-—-- = hi -— + --~--- + ---- XI-8 
ax k dx x (s+x) 
dk 
The term -— may be written as: 
ax 


dk dk do... 


ax ao dx 
and with the aid of Equation XI-4, there results: 
dk 86-98 dk 
2 ae es X19 
Substitution from Equations XI-4, XI-§ and XI-9 into 


Equation XI-e yields: 
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x 
i +m- 1] s(nea)(--2--)} XI-10 


B+ & 

This is a startling result. If the surface is iso- 
thermal and lacks a hydrodynamic starting length, the sub- 
stitution of previously aevsrmined values for m and n with 


turbulent boundary layers, yields: 


h 
-- = 0.80 Siam 
Dn 


This expression shows that the local coefficient at a point 
x is eighty per cent of the value of the mean coefficient 
to the point, 

A comparison with Latzko's theoretical results for the 
flat plate and turbulent boundary layers may be made. 


Latzko's result for the local surface coefficient is: 


2 
| 


Oe 

= 0.0265 Seve (#2 | XI-le 

Or OLS aL 
O 

= hd, 


Latzko's result for the mean surface coefficient is: 


J Oa 
q = 0.0356 8 cV,0,x ca XI-13 
Ve ox 
= DyPn* 
The ratio of Equation XI-le to Equation XI-13 yields: 
h 0.0285 
Hes 0.0356 


which is the same as Equation XI-1ll. 
Equation XI-10 may be put into a more useful form. For 


a small change in the fluid temperature, it may be assumed 
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that: 


{ 
x + BY (ty) 


( 
where ts comes from Equation VIII-22 as: 


k=a 


‘ 
ter = te + one 


so that: 
k =ea, + b(t +0 /2) XI-15 
Differentiation of Equation XI-15 yields: 
es 
de, a 


Reference to the NBS-NACA Tables (27) permits evalua- 


tion of b,. The data may be given as: 


Temperature, °F k, B nyt g¢7t pol 
62 1.472 x 1072 
98 1.562 x 107 

5 ee 


Tnese data give b; to be 2.5 x 10 “B Re F 

The minimum value of k in the experiments is not less 
than that at 62°F. Examination of the data shows that the 
maximum value of 6 - 0, may be 5 F, approximately. With 


6-6 dk 
these extremal values, evaluation of the term ------ a 


de, 
yields 0.004 as a maximum. Therefore, the term may be neg- 
lected when compared with the developed value for h/h,, in 
Equation XI-ll. 
Insertion of other known values into Equation XI-10 
yields: ) 
h/h, = 1 - @,/0(0.09 + 0.11x/L) XI-16 


It is of interest to test Equation XI-16 with 
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experimental data. The comparison between calculated val- 
ues and experimental results may be made in Table XXV. Ex- 
amination of the table shows that the experimental result 
ls greater than the calculated value in almost all cases, 
The maximum difference between the experimental value and 
the calculated value is about 12.9 per cent, based on the 
experimental value. Thus Equation XI-16 would give conser- 
vative results for h. 

The agreement between experimental and calculated re- 
sults is very good for values of s/L smaller than 0.95. 

The majority of differences between experimental and cal- 
culated results are of the order of five percent of the 
experimental result. 

When consideration is given to the nature of the app- 
roximations required for the calculation of the surface 
coefficients from the basic data, the agreement between 
experimental and calculated results is considered to be 
satisfactory. 

A surface configuration function for the local surface 
coefficient may be obtained. If it is assumed that the 


surface is isothermal, then Equation XI-16 yields. 


t 
(Ny) s = (Ny), (0,91 =O. 2677 1h) 
Substitution for (Ni) from Equation X-2 yields: 
' : 0.91 1.91 
(iy), = 0.0307 (tp)2°® [0.91¢x/1)°°% - 0,22(¢x/2)?* 9] 


XI-17/ 
This equation may be put into a form wherein the term in 


the brackets becomes unity when x is made equal to L, if 
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the constants in the brackets are multiplied and the equa- 
tion constant is divided by 1.25. The result is: 


(Nyx = 00246 (Npe)P?® [1.138(x/1)°* 9 0.136(x/L)2*92] 


Xue 
ad, I-18 


0.0246 (Nre)p [2..238(x/L-0°97_ 0.138(x/L) °° 9+] 


eran: 


A comparison of surface configuration functions is 


(Nya) 


made in Tables XXVI and XXVII. Equation II-& in combina- 
tion with Equation III-2 remains unchanged if c.. is divided 
and § is multiplied by 1.25. In this way the configuration 
function of Jakob and Dow may be made comparable to the 
present configuration function. The method has been shown 
in passing from Equation XI-17 to Equation XI-18. This pro- 
cedure has the effect of permitting the function to assume 
the value of unity when x becomes equal to L. 

Examination of the tables in which extrapolated values 
are given in parenthesis shows that agreement among the 
various surface configuration functions is satisfactory for 
values of x/L greater than 0.3. Greater discrepancies occur 
for the smaller values of x/L. Except for x/L equal to zero, 
the functions of Rubesin, Jakob and Dow, and the present work 
yield acceptable results for small values of x/L. Obviously, 
extrapolation of Jakob and Dow's empirical form to very small 
values of x/L is not possible because this form does not con- 
sider that the function will become infinity as x/L app- 


roaches zero, 
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TABLE XXV 
COMPARISON OF LOCAL SURFACE COEFFICIENTS 
Data selected for (Noe dy equal to or greater than 5 x 10?. 


Series A Runs 


mem x, in. h/h,, Difference, per cent of 
No Calc. Exp » experimental value 
13 16 0.825 0.879 6.2 
20 0.820 0.907 9.2 
ol 0.813 One 12.9 
14 12 0.829 0.869 4.6 
16 0.825 0.885 6.8 
20 0.822 0.903 9.0 
ey 0.815 0.916 110 
15 12 0.832 0.851 2.2 
16 0.828 0.83 6.1 
20 0.822 0.897 8.3 
ey 0.81% 0.929 12.5 
Us S 0.835 0. S45 dae 
12 0.833 0.859 3.0 
16 0.828 0.619 5.4 
20 0.823 0.89 .0 
ey 0.615 0.914 10.8 
Series C Runs 
& 12 0.871 0.86 ~0.7 
16 0.859 0.868 2.8 
9 & 0.881 0.889 153.0 
12 0.871 0.889 220 
16 0.859 0.902 4.8 
10 S 0. &&0 0.916 3.9 
12 Mo 93 0.875 0.5 
16 0.60 0.899 4.3 
td yy 0.890 0.938 5.1 
S 0.881 0.881 O 
12 0.870 0.896 2.9 
16 0.659 0.901 4.7 








Run 
No. 


Pc 


13 


14 


15 


10 


ie 


13 


xX») 
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TABLE XXV, Continued 


Series C Runs, continued 


Calc. 


0.891 
0.880 
0.8/0 
0.859 


hy 


69) 
@ 
ae 
}- 
@ 
1¢8) 


Exp. 


0.934 
01364 
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TABLE XXVI 
COMPARISON OF SURFACE CONFIGURATION FUNCTIONS 


For Local Surface Coefficients and Turbulent Boundary Layers 
Extrapolated values in Parentheses 


(x/L)  Rubesin Maisel and Jakob and Present (s/L) 
Work 


Equ. III-5 Sherwood Dow or 
Equ. III-%3 Equ. III-e Equ. XI-19 

0.0 mars ( of ) Leiae) Cae, 1.0 
Or 2 1.518 1.318 bee 1. 383 0.9 
2 a ci ger? bey 1.281 0.8 

. : 1.200 1.22 O. 
074 1.183 1.128 1.130 Rely A: 
0.5 1.136 1.09% 1.090 1.138 7 
017s tloeg, «(120573 i083 slow? 
0.8 1.042 Hees 1.042 1.049 0:2 
0.9 eee 1. tant 1.025 1.023 Onn 
ir, 1.000 (1.000 (1.000) 1.000 0.0 


TABLE XXVII 
DEVIATION OF SURFACE CONFIGURATION FUNCTIONS 
For Local Surface Coefficients and Turbulent Boundary Layers 
Body of table is the deviation in per cent, based on Rubesin's 
Extrapolated <p Parentheses 


(x/L) Present Jakob and Maisel and (s8/L) 


Work Dow Sherwood 

0.0 C2) (100) [) Ie 
O.1 <- 8.9 = J), -~ 13.2 0.9 
0.2 —- 4.5 a eee) 9.0 0.8 
Ss =n) 5 o.7 
ol = als a as = 422 O 

0.5 + 0,2 al - 3,4 a 
O. + 0.4 = 2as (~ 2,2) oO. 

0.7 +0.8 -~ 1.5 (- 1.2 0.3 
0.8 +0.5 O : Oa7 0.2 
Oo + 0.2 + 0.3 = O23 O.1 
1.0 O ( 0) (O 0.0 








CHAPTER XII 
THE EFFECT OF SURFACE CURVATURE 

The experimental results with respect to curvature are 
conveniently compared by aSsuming unit Reynolds number, 
zero hydrodynamic starting length, and unit total length. 
This procedure serves to make Nusselt number equal to oo of 
Table IV and the influence of curvature is reflected there- 
in, 

The available experimental results for curved surfaces 


with air at atmospheric pressures and temperatures are: 


t 2 

Source Cy Kyi - 
Jakob and Dow _ 0.0280 18.46 
Present Work 0.0307 38.46 


Jakob and Dow have given a correction factor, a result 
of theoretical reasoning, which is a multiplier for Latzko's 
result for the flat plate. The factor is: 

F = 1 + (3/10) bK. XII-~1 
It is noted that F is a linear function of the curvature. 
Therefore, it is assumed that C. may be given in the form: 
c. 


o = conn byK . XII-e 


Substitution of the two experimental results enables solu- 
tion of Equation XII-2. The solution yields: 
t = 
Co = 0.0255 + 1.35 x 10 : K , XIIT-%3 
4 
wherein C, is 0.0255 and by is 1.35 x 10 ft. 
Extrapolation of Equation XII-3 to zero curvature 
t 
yields C, as equal to 0.0255. It may be recalled that 


Latzko's value for the flat plate is 0.0253. These results 


14s 
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are illustrated in Figure 38. 
The theoretical reasoning of Jakob and Dow may be sun- 


marized as: 


w= 2 1r TV.c @,b 6 (eaee= a pean (1 + * (meee 
(n+1)(2n+1) r (nt2)( 2n+2) 
XII-4 
This equation comes from considering a heat and mass balance 
for an elemental slice through the boundary layer in which 
it is assumed that the velocity distribution is given by: 
v = (y/b) vy > XII-5 
and that the temperature distribution is similar. 
Jakob (30) has shown that the exponent n in Equation 
XII-5 may vary from 1/6 to less than 1/10 for flow inside 


tubes. It is of interest to determine how similar changes 


may affect Equation XII-4. The computation yields: 


n (n+1)(2n4+1) 

(nt1)(2n+1) (n+2)(2n42) 
1/6 Oeaaloyyal O. 308 
aL 0.0973 O. 200 
1/8 0.0889 0.294 
1/9 0.0818 0.289 
i 0 0.0758 0.286 


This computation shows that a larger exponent yields greater 
heat transfer. 

It is noted from Table IV that Colburn's correlation 
gives the smallest experimental result for the flat plate, 
in which the curvature is zero. Latzko's result is based 
upon the exponent n being 1//7, as suggested by von Karman 


(31). Comparison of these two results shows that Colburn's 
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value is 1.265 times greater than Latzko's value. Such a 
result requires the exponent n to be about 1/4. 

On the other hand, the effect of the exponent on the 
factor F is seen to be small in comparison with the effect 
of the exponent on the flat plate result. The Jakob and 
Dow result for a cylinder 1.% inches in diameter is about 
ten per cent higher than Latzko. If the curvature has no 
effect, then this increase is of the same order of magni- 
tude given by n equal to 1/6. 

The present result is about twenty-one per cent higher 
than that of Latzko. If the curvature has no effect, then 
n should be about 1/4 as with Colburn. Thus the results 
for the flat plate are not compatible with those for curved 
surfaces for any reasonable change in the value of the ex- 
ponent, 

Equation XII-4 shows that an increase in the boundary 
layer thickness causes an increase in the heat transfer. 

It is possible that the inter-relation between the boundary 
layer thickness and the exponent n is such that the net 
effect on the change in heat transfer between the plate and 
a cylinder is zero. However, by the Jakob and Dow result, 
it is seen that the inter-relation reaches a minimum value 
somewhere between the present result and that for the plate 
as given by Colburn. Such a condition may not occur for 
reasonable changes in the exponent because it has been shown 
that the exponent would have to change from 1/4 to 1/6 and 


back to 1/4 with corresponding changes in the boundary layer 





152 


thickness, 

AS a consequence of this discussion, it is concluded 
that the curvature is a variable of great importance. 
Reasonable changes in the value of the exponent or of the 
boundary layer thickness do not account for the experimental 
results with cylinders. 

The influence of curvature may be added to previously 
determined results, yielding: 


= (0.0255 + Dy/T)(Npeoo(x/L)o* + 


(Nw) x Gi 


and 
wy) p, = (000255 + dy/t) (Nag) pO x/L) O°? x1 


in which b, has the value 1.35 x 107* ft and r is 1/K ft. 


(N 





CHAPTER XIII 
SUMMARY AND CONCLUSIONS 

1. A heat transfer element was constructed having a 
nominal diameter of 0.624 inches. Cylindrical nosepieces 
of various lengths having hemispherical nosetips were att- 
ached to the heat transfer element. The nosepiece lengths 
were varied in the proportion 55:1, and the thermal lengths 
were varied in the propertion 7:1, so that the ratio of 
thermal length to total length varied from 0.141 to 0.983. 

2. The heat transfer element was placed in an air 
stream so that the surface generators were parallel to the 
flow. The stream velocity was varied from 9.9 ft sec tto 
ae.S ft sec corresponding to a variation in Reynolds 
number from 2.6 x 10° to 2.4 x 10° based on a defined total 
length as the characteristic length. Transition to turbu- 
lence started at Reynolds number of 50,000 and was completed 
for Reynolds number as low as 3 x 10? and as high as 6 x 10. 

3. The result of measurements was obtained for turbu- 
lent boundary layers only. The result could be expressed 
as: 


(Nyy) g = 0+0307 (Npg)p °(x/L)°°7* , oe 


Nu) x 


where (N is based on the mean surface coefficient of heat 


Nu? x 
transfer and the thermal length x as characteristic length, 
and the Reynolds number based on a defined total length L 

as characteristic length. The above equation was found to 


be valid for variable thermal length as well &s variable 


hydrodynamic starting length. 
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4, An equation for the Nusselt number based on the 
local surface coefficient of heat transfer was developed 
from Equation X-e by differentiation. The developed equa- 
tion was tested with experimental data. The developed 
equation gave results lower than the experimental measure-— 
ments varying from zero to le.9 per cent, with an average 
variation of five per cent. 

5. The total length of the heat transfer element and 
attached nosepieces was such that an unsupported overhang 
of three feet occurred. Vibration and cross flow studies 
were made and it was concluded that the experimental results 
were not influenced by these factors, 

6. The results of a previous investigation by Jakob 
and Dow were used in conjunction with the present results 
to determine the influence of curvature. The final result 
could be expressed as: 

(Ny), = (0.0255 + d,/r)(N, Pe O(x/L)OF2 , XII-6 
by 


where by, = 1.35 x 10°" ft and r is the radius of a section 
of the surface perpendicular to the air stream. 

7. Equation XII-6 coupled with the theory of Latzko 
for the flat plate leads to the conclusion that the exper- 
imental results for the flat plate reported by different 


previous sources are in error and that Latzko's result is 


more nearly correct for the flat plate. 














APPENDIX I 
VOLT BOX MODIFICATION 

Modification of a Leeds and Northrup Co. volt box, No. 
418271, was undertaken so that all direct current electrical 
potential differences likely to be encountered could be 
measured with a potentiometer. The portable precision po- 
tentiometer at hand was capable of measuring potential diff- 
erences up to eighty millivolts, whereas the volt box was 
designed for use with a 1.5 volt potentiometer, 

A volt box consists of a high resistance coil provided 
with taps to accomodate various impressed potential differ- 
ences. An additional tap was provided so that the precision 
potentiometer could be used. 

The volt box was provided with terminals for impressed 
voltages of 15, 150, and 300 volts. These terminals are 
designated as Eis: E150: and E300: respectively. The init- 
jal measuring terminal is designated as o5 and the addit-— 
ional measuring terminal is designated as ©9075" 

After modification, the volt box was tested by measur- 
ing the various resistances with a Wheatstone bridge. These 


measurements were: 


Coil Portion Resistance, ohms 
Ro 075 > 15.09 
45 300.2 
Ris 4000 
R150 30 ,Q00 
R200 60,000 
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where the subscripts correspond to the voltage terminals. 

An additional test was made by connecting 1.5 volt dry 
cells in series to the volt box. Simultaneous observations 
were made at the measuring terminals. 

The test data are given in Table I-l. These data show 
that the ratio €1 5/€0.075 may be 19.92. The probable error 
may be obtained from the formula: 

re = 0.8453 [n°(n - 1) | -2/? pa 
where r, is the probable error of the ackn, 
n is the number of independent observations, 
x. is the difference between the i*® observed value 
and the arithmetic mean, 

The comoutation yields: 

te AZ 0.6453 (256 x 15) 71/2, 0.247 
== 0.0034 


The test data may be converted into a more useful form 


as: 
E300 = Sates = 398424 O76 
£150 = 100e) 5 = LIFECD 075 
Eis = 106) 5 = 199-28 075 


where the e are the potential differences at the measuring 
terminals in volts. 


These results may be compared with that obtained with 





litandbook of Engineering Fundamentals. Ovid W. Eshbach, 
Editor. Wiley Engineering Handbook Series. First ed. Elev- 
enth printing. John Wiley and Sons, Inc. New York. 1947. 
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the Wheatstone bridge as; 
©) 5/°0.075 = 300.2/15.09 
= 19769 . 
It may be noted that the maximum current that flows in 
the volt box at rated impressed potential differences is: 
I = 300/60,000 = 150/30,000 = 15/3000 
= 0.005 amperes. 


TABLE I-1 
VOLT BOX TEST DATA 


Number of Terminal a e Ratio Dithi. 
Dry Cells 1.5 0-0/5 from 
in series mv mv €1 5/£0.075 mean 
1 300 7.923 0.398 19.907 0.017 
1 300 7.925 0.398 19.912 0.012 

2 300 15.740 0.790 19.924 0.0 

2 300 15.740 0.790 19.924 CROs 
3 300 23.53 1.182 19.90 0. Oli 

; 300 23.53 Pee 19.92 0.0 
300 gu. 5D 1.575 1 OD 0.019 
yy 300 1.35 1.575 19.905 0.019 

4 150 ae 3.149 19.911 0.01 

4 150 62.69 3.149 19.908 0.01 

3 150 47.06 2. 363 19.915 0.00 

3 150 47.0 2.363 19.920 0.00 
2 150 31.4 1.578 19.937 0.013 

2 150 31.46 we5 79 19.924 0.0 
1 150 15.840 0.792 20.000 0.076 
1 150 15.845 0.792 19.956 0.032 
TOtads...+. 318.779 0.247 


Arithmetic mean.. 19.924 





APPENDIX II 
CALIBRATION OF WESTON SHUNT 

A Weston Electrical Instrument Co. shunt, number 33956, 
was used to measure the flow of electric current where nec- 
essary. This shunt was provided with five-ampere, twenty- 
ampere, and fifty-ampere terminals. Of these, the five-ampere 
and the twenty-ampere portions were tested, 

The method of testing used was to connect this shunt in 
series with @ standard resistance and then allow an electric 
current to flow through them. The potential drop caused by 
each shunt was measured with a potentiometer at various flows 
of current. The source of current was a storage battery. 

The standard resistance was a Siemens, number 2006346, 
rated for fifty amperes in air, made of Manganin, and having 
a resistance of 0.01 ohms at eO deg. C. 

For the five ampere range of the Weston shunt; using the 


formula: 


r 7 0. S45 3x [n?(n-1)| Sx, 


where Ip is the probable error of the mean, 


n is the number of independent observations, 


,th 


X is the difference between the observed value and 


i 
the arithmetic mean, 


' there results: 1/2 
rs mz 0. S453 x [1600( 39) | x 0.0097 


a ee 


luandbook of Engineering Fundamentals. Ovid W. Eshbach, 
Editor. Wiley Engineering Handbook Series, first edition 
Oe printing. John Wiley and Sons, Inc. New York. 1947. 
pg e~leb. 
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oper sie =—= 0.00003 
The value of the ratio used is 0.4340. 
For the twenty ampere portion, the value of the ratio 
used is 1.844, 
Since the electric resistance of the standard shunt is 
O.O01 ohms, using Ohm's Law, it follows that: 
for the five ampere portion of the Weston shunt, 
I = 0.04340 x EO» and 
for the twenty ampere portion of the Weston shunt, 
I = 0.1844 x E 
where I is the electric current in amperes, and 
Ew is the potential difference across the Weston shunt 


in millivolts. 








E 1 


1 


Totals 


w» if 
millivolts 
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TABLE II-1 


CALIBRATION OF WESTON SHUNT 


Eg, i Ratio 
millivolts E,/Ew 
Five ampere portion 
2.099 Oa 
2.09 0.433 
2.09 0.4335 
2.095 0. 1339 
2.97 0.433 
2.97 0.4338 
omen 0.4339 
Oe 0.4339 
6.018 0.4342 
6.020 0.4340 
6.021 0.4340 
6.022 0.4340 
9.23 0.4331 
9+25 oe 
9 te 0.4342 
12.07 0.4342 
12.07 0.4340 
12.07 0.4340 
15.13 0.4339 
15.12 0.43 
15.1 0.4345 
15.19 0+ 990 
18.29 0.43 
18, 30 0.4342 
18. 30 0.4347 
al. Se 0.4347 
21. 36 O.4 
21. 38 0.4340 
21. 38 0.4340 
25.99 0.4343 
25.98 0.4342 
30. 80 0.434) 
30,82 0.4342 
30.81 0.4341 
30.81 0.434 
733599 


DiI. 


0.00012 
0.00968 
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TABLE II-1, Continued 
CALIBRATION OF WESTON SHUNT 


Ey» mv Eg» mv Ratio, E,/E, 
Twenty ampere portion 
16.96 sey 1.344 
16.96 31... 1. S44 
a Ui Se a 20.54 1.54 
ries 20.52 1.84 
5.9/1 TiOve 1.344 
5.97e 11.014 1.344 
Room temperature: (For both portions) 
Ab start... ..... a np 


Diy WH. we ee eo - (7 





APPENDIX III 
CALIBRATION OF THERMOCOUPLES 

The seven thermocouples were calibrated simultaneously. 
The hot junctions were in an atmosphere of condensing steam 
in a hypsometer. The cold junctions were at the temperature 
of a mixture of melting ice and water at atmospheric press-— 
ure. The hypsometer, vented to the atmosphere, was supplied 
with the steam from distilled water and was equipped with a 
water manometer for indicating interior pressure excess. 


The data and results are as follows: 


Thermocouple emf, millivolts Equivalent temp., °F. 
No. 
1 4.258 ran i 
2 4.258 Fadl ied Oya 
3 4,258 elec 
Four in series 17.00 adele © 


The computed temperature of condensing steam for a corrected 
barometric reading of 743.1 mm of mercury is 210.9 °F, using 


the Keenan and Keyes formula (25) 
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APPENDIX IV 
CALROD HEATER SURFACE TEMPERATURE EXPLORATION 

The determination of the uniformity of heat generation 
along the length of the Calrod heater was most important. 
Direct measurements of this characteristic could not be 
made, 

A convenient indirect method is provided by measurement 
of surface temperature. If a heated element, surrounded by 
still air at a uniform temperature, is supported horizon-— 
tally, then heat is transferred by free convection and radia- 
tion. The surface assumes an equilibrium temperature when 
the steady state is reached, 

If the heat generation is uniform, the surface temper- 
ature will be the same over the whole length. Likewise, the 
surface coefficient of heat transfer remains constant, 
Variations in surface temperature along the length may be 
expected if the heat generation is non-uniform. 

Study of the theory of heat transfer by free convection 
shows that the error made by assuming that the coefficient 
of heat transfer remains constant for small surface temper- 
ature variations is negligible. Also, if the surface temp- 
erature is higher at one place than at another place on the 
element, then the heat transferred from the one place is 
greater than that from the other place. It follows, there- 
fore, that the heat released from the coil is greater at the 
one place than at the other. Thus the measurement of the 


surface temperature becomes a means to determine the 
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uniformity of heat generation. 

A knowledge of the specific surface temperature is not 
required. The temperature distribution along the element 
is sufficient. 

The measurement of the surface temperature was accomp- 
lished by suspending a movable thermocouple from the heater. 
This procedure causes an interference with the heat transfer 
process at the point of application of the thermocouple, but 
it is reasonable to assume that the influence of the inter~- 
ference is the same at all points of application. 

The thermocouple, non-standard, was made of: 

Chromel P, No. 40 Brown and Sharp Gauge, 0.0031 inches 
Giameter, 

Alumel, No. 40 Brown and Sharp Gauge, 0.0031 inches 
diameter. : 


The calibration data for a thermocouple made of these 


wires are; 


Temperature, Millivolts Error, @F 
oF No. 40 ga. No. 18 ga. No. ga. 
i 0.998 0.998 O 
161 Coo] fe A at E 
227 sta 7 4, 40 + 
668 14.450 14.430 +] 


The No. 18 gauge thermocouple follows the Leeds and Northrup 


Co. tables, Standard 31031. 


seers evenmetnven tn ear at Ee PSEC TOSSES SL SS NN 


lpr. M. Spielman, Associate Research Engineer, Department of 
Mechanical Engineering, Illinois Institute of Technology, 
was kind enough to furnish these data and to show me the 
method he evolved for butt-welding these wires. 
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In service, the hot junction was located so that the 
welded joint was at the top surface of the Calrod heater 
being in contact with the surface for a distance of one- 
half the circumference of the heater, 

A block of wood suspended eighteen inches below the 
heater served as the site for the cold junction for this 
System, The chromel wire was soldered to a copper—constan- 
tan thermocouple and the alumel wire was soldered to a cop= 
per wire. These wires were fastened and glued to the wood 
block. 

The temperature of melting ice in water was the cold 
junction reference temperature for the copper-constantan 
thermocouple. 

Two series of runs were made with this thermocouple 
System. One series was made with an electrical energy in- 
put of 125 watts and the other at 41.2 watts. 

The recommended method for compensating for the cold 
junction temperature when the cold junction is at a temp- 
erature different from standard is to convert the observed 
temperature at the cold junction to the millivolt equivalent 
of the thermocouple? In the present case, the cold junction 
is at a higher temperature than the reference standard, 
hence the millivolt equivalent is added to the observed po- 


tential difference to obtain an argument for use of the 





lapplications of Chemical Engineering. H. McCormack, Editor. 
first edition. Second printing. D. Van Nostrand Co., Ine. 
New York. 1940. 
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conversion tables. 

It may be recalled that the cold junction temperature 
of the chromel-alumel thermocouple was measured with a 
copper-constantan thermocouple. To reduce the commission 
of errors of conversion, the data from the Leeds and North- 


rup Co. tables, Standard 31031 may be: 


Temperature, Millivolts, 32 °F reference 
a Copper-constantan Chromel—-alumel 
100 1.52 aoe 
she, 1.29 1.29 
&O 1.06 1.06 
{0 0.83 0. 54 


Inspection of these data shows that the two kinds of thermo-— 
couples may be considered to be thermometrically equivalent 
for the range of temperature indicated. 

The temperature at the wood block was considered to be 
the temperature of the air passing across the heater surface, 

The observed data are given in Tables IV-1 and IV-e. 

It is seen that one end of the heater was hotter and that 
the heat release from the coil was non-uniforn. 

A reasonable estimate of the temperature distribution 
along the heater may be obtained by assuming that the temp- 
erature variation is a linear function of the length. This 
estimate may be obtained by computing a least squares line 
to fit the data with the middle twenty-four inches as a 
basis. 

Computation yields the following: 

ee 367.1 - 0.43 2, F, and (Part A) 


a) 160.4 - 0.19 z, PF (Part B) 


B 








169 
TABLE IV~1 
CALROD HEATER SURFACE TEMPERATURE 
Part A: 125 watts input 


Distance along Surface temp. minus Variation from 
heater, in. room temp., F mean, F 
Run 1 Run 2 Run 1 Run 2 
7 316 316 - S820 -— eT 
: ; - : 3 7 Se 
10 36 366 + 12.0 + 6,3 
11 365 367 *#11.0 + 7.3 
Ve 362 359 + 6&0 <= 0.7 
i 356 356 + 2.0 = 3.7 
v 356 356 “<0 = Fa 
1 361 361 + 7.0 + Eve 
ya 363 365 + 9.0 + 5,3 
1g 362 369 + 60 + 9,3 
19 362 366 + 8.0 + 6,3 
20 357 368 + 3.0 + 8,3 
21 360 362 + 6.0 + 2,3 
22 357 364 + 3.0 + 4,3 
2 351 360 -~ 3.0 + 0.3 
2 352 361 - 2.0 + 1.3 
e 39 DoD = a - 4.7 
2 il = o e 
f! 326 360 + 2,0 + 0,3 
2g 351 358 _ R°0 - 1.7 
30 350 39 - 4.0 = 0.7 
31 351 301 =~ 50 + ae 
32 351 360 - 2.0 + G3 
3 7 56 - (.0 (= Wee 
3 345 3 9g as 9.0 — QO. / 
37 341 3 - 13.0 -15.7 
38 Bib. 3 ~ 10.0 - 15.7 
9 346 346 ~ 8&0 -13.7 
41 311 Bpial ~ 43,0 - 48,7 


Arithmetic mean of 
middle twenty-four 
6). eee 354.0 359.7 





1/0 
TABLE IV=2 
CALROD HEATER SURFACE TEMPERATURE 
Part B: 41.2 watts input 


Distance along Surface temp. minus Variation from 
heater, in. room temp., F mean, F 
Run 1 Run 2 Run 1 Run 2 
1 132 140 - 25.0 - 24,7 
& 150 Bi ~ 7.0 - 6.7 
157 156 0 + ‘leas 
10 160 160 op he: + Fo 
11 161 157 + 4,0 + 2,3 
le 159 15 + 2,0 + 0.3 
al 15 15 0 =- oy 
1 15 154 hae - QO./7 
1 156 152 + 120 - 2.7 
1 157 154 0 - O./7 
17 158 158 ae. O + ee 
18 159 159 + 2.0 + ae 
19 161 159 + 4,0 + 4,3 
20 162 158 + 20 + 3,3 
21 161 158 + 4,0 + 3.3 
ee 159 156 + Gano + i233 
2 159 155 + 2,0 + O26 
2 154 155 - 3.0 + 0.3 
27 157 155 O + 0.3 
28 15 155 - 2,0 + 0.3 
29 15 15 - 1.0 + 0.3 
30 155 L5 ~- 2,0 =- O37 
31 155 154 - 2.0 - 0.7 
3 2 IDS 8, - 3.0 - “Lae 
3 15 pS at - 2,0 - 8.7 
3 153 151 ~ 0 “ Agll 
G2 1 14 - 10.0 - 10.7 
Ky 13 131 - 23,0 - 23.7 


Arithmetic mean of 
middle twenty-four 


inches 


_ ee 2 yO 154.7 
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